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Properties of GRB Lightcurves from Magnetic Reconnection 

Paz Beniamini^ and Jonathan Granot^ 

ABSTRACT 

The energy dissipation mechanism within Gamma-Ray Burst (GRB) outflows, 
driving their extremely luminous prompt 7-ray emission is still uncertain. The lead¬ 
ing candidates are internal shocks and magnetic reconnection. While the emission 
from internal shocks has been extensively studied, that from reconnection still has few 
quantitative predictions. We study the expected prompt-GRB emission from magnetic 
reconnection and compare its temporal and spectral properties to observations. The 
main difference from internal shocks is that for reconnection one expects relativistic 
bulk motions with Lorentz factors T' > a few in the jet’s bulk frame. We consider 
such motions of the emitting material in two anti-parallel directions (e.g. of the recon¬ 
necting magnetic-field lines) within an ultra-relativistic (with T 3> 1) thin spherical 
reconnection layer. The emission’s relativistic beaming in the jet’s frame greatly af¬ 
fects the light-curves. For emission at radii Rq < R < Ro + Ai? (with T = const) the 
observed pulse width is AT ~ (i?o/2cr^) max(l/r', AR/Rq), i.e. up to ~ T' times 
shorter than for isotropic emission in the jet’s frame. We consider two possible mag¬ 
netic reconnection modes: a quasi steady-state with continuous plasma flow into and 
out of the reconnection layer, and sporadic reconnection in relativistic turbulence that 
produces relativistic plasmoids. Both of these modes can account for many observed 
prompt-GRB properties: variability, pulse asymmetry, the very rapid declines at their 
end and pulse evolutions that are either hard to soft (for T' < 2) or intensity tracking 
(for r' > 2). However, only the relativistic turbulence mode can naturally account 
also for the following correlations: luminosity-variability, peak luminosity - peak fre¬ 
quency and pulse width energy dependence / spectral lags. 


1. Introduction 

Gamma Ray Bursts (GRBs) are the most extreme explosions in nature. They are powered by 
ultra-relativistic jets (with Lorentz factors T > 100) and have huge isotropic equivalent luminosi¬ 
ties (Liso ~ 10^^ erg/s). A good part of the outflow energy is thought to be dissipated at large 
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distances (R ~ 10^^ — 10^^ cm) from the central engine, where it efficiently radiates, predomi¬ 
nantly in 7 -rays that escape the emission region and eventually reach us. One of the major open 
questions in GRB research concerns the outflow composition. In particular, whether the energy 
is carried out from the central source to the emission region predominantly as kinetic energy - a 
baryonic jet (Shemi & Piran 1990), or as Poynting flux - a highly magnetized (or Poynting flux 
dominated) jet ( Usov [1992 ; Thompson |1994[|Meszaros et al. |1997t|Lyutikov & Blandford [2003 


Granot et al. 2015). Recently, Beniamini & Piran (2014) have shown that the question of outflow 
composition is essentially tied to another major open question in GRB research - the nature of the 
prompt emission mechanism. In particular, one zone emission models in magnetic jets naturally 
produce efficient synchrotron radiation and predict fluxes in the X-ray and optical bands in excess 
of the observed ones during the prompt. However, the magnetization might still be high up until 


the emission region, and then sharply drop within the emission region itself (Drenkhahn & Spruit 
2002 [ Lyutikov & Blandford |2003] Sironi et al. |2015 ), if the dissipation that causes the emission 
is driven by efficient magnetic reconnection. 

Magnetic jets can be powered by the rotational energy (ultimately arising from gravitational 
energy) of a magnetized torus around the black hole ( Meszaros et al. ||1997 ) or of a highly mag¬ 
netized millisecond neutron star ( Usov ||1992 ; Spruit 1999| ). Initially highly-magnetized jets are 
favored on energetic grounds, as modeling of GRB central engines that rely on accretion disks 
suggest that their power is significantly larger than that of thermally driven outflows powered by 
neutrino-anti neutrino annihilation (see e.g. Kawanaka, Piran, & Krol ik 2013[ ). Moreover, a 
high value of the magnetization parameter a (the magnetic to particle enthalpy density or energy 
flux ratio) near the source can help prevent excessive baryon loading that may prevent the jet from 
reaching sufficiently high Lorentz factors far from the source, at the emission region. Such initially 
highly-magnetized jets may naturally lead to magnetic reconnection. In a striped wind magnetic 
field configuration ( Coroniti|1990 ), with either periodic or stochastic flipping of the magnetic field 
direction near the source, reconnection at large distances from the source has a natural preferred 
direction. Moreover, for large cr-values just before the dissipation region, reconnection will lead 
to local relativistic bulk motion of the outgoing particles away from the reconnection sites (in the 
bulk frame of the jet). Therefore, reconnection models may be able to overcome the difficulties 
present in an isotropic arrangement of the magnetic fields regarding overproduction of optical and 
X-ray radiation (as compared with observations) due to synchrotron emission. 

As mentioned above, synchrotron is likely the dominant source of emission in magnetic jets. 
Many studies have considered the implications of synchrotron radiation being the source of the 
prompt emission in GRBs ( Katz|1994 Rees and Meszarosj 1994 ; [Sari et al.||1996||199^ ; Kumar & 
McMahon 2008^ Daigne|2011 Beniamini & Piran 2013). One of the greatest challenges for recon¬ 
ciling this model with the observations is related to the low energy spectral slope. The synchrotron 
fast cooling (which is the likely cooling regime in order to achieve the required large radiative 
efficiencies) photon index below the peak, dN/dE oc i?"® with as = —1.5, is inconsistent with 
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the average observed slope of = —1. Although it may be possible to aehieve softer speetra, 
as < —1.5 (for instanee by having the peak frequeney fiuetuate on a shorter time-seale than that 
for whieh the speetra are obtained, and in effeet smearing the transition from below to above the 
synehrotron peak), it is not at all trivial to inerease in these models. However, 90% of all GRB 
have a > —1.5 ( [Preeee et al. ||2000HGhirlanda et al.|2002[|Kaneko et al. |2006[|Nava et al.|2011| ). 
In about 40% of the GRBs a > —2/3 (Nava et al. 20111, whieh is impossible for synehrotron, even 


in the ease of slow eooling. This problem is known as the synehrotron “line of death” ([Preeee et al. 


19981. However, a refined time-dependent speetral analysis in several bright bursts (Guiriee [2012 ) 


has suggested that GRB speetra may be better fitted with a multi-eomponent model instead of the 
elassieal “Band funetion”. In this model, the Band funetion eomponent dominates the total energy 
but is aeeompanied by weaker power-law and blaek body eomponents. The lower energy speetral 
slope of the Band funetion is softer than in Band only fits, and is eonsistent with slow eooling 
synehrotron (i.e. the synehrotron “slow eooling line of death” problem is removed). Reeonneetion 
models may naturally allow for eontinuous heating of the eleetrons while they are emitting, in this 
way produeing low-energy speetral slopes eonsistent with slow eooling while maintaining large ra¬ 
diative effieieneies ( jGhisellini and Celotti||1999[ [Kumar & MeMah^|2008[ |Fan||2010[ jBeniamim 
fc Piran 12014] ). 


An anisotropie emission model has been suggested to operate in the afterglow phase of GRBs 
(Beloborodov et al. 2011| ). These authors argued that synehrotron radiation is almost eertainly 
anisotropie, and suggested that this anisotropy eould help explain highly variable signals observed 
in some X-ray afterglows (sueh as the rapid deeay phase or X-ray flares). In this work, we explore 
an analytie model in whieh the emitting eleetrons are anisotropie in the jet’s bulk frame represent¬ 
ing the bulk motion of the jet, or more speeifieally, are isotropie in their loeal eenter of momentum 
frame, whieh moves relativistieally relative to the jet’s bulk frame. This eould be applieable to 
reeonneetion from a striped wind where eaeh pulse in the light-eurve of the prompt emission is 
due to some maeroseopie reeonneetion event ( Kumar & Crumle}^|2015| ). Moreover, it ean also 
be applied in a broader eontext of anisotropie emission, sueh as had been suggested in the “jet in 
jet” model ( [Levinson et al. |[T993[ [Pedersen et al. [[T998| [Frail et ar~[[2000[ ) and the relativistie 
turbulenee models ( jLyutikov & Blandford [2003] [Kumar & Narayan |2009[ [Lazar et al. |2009[ ). 


The paper is organized as follows. In § [^ we present two different physieal seenarios that 
ean allow for anisotropie emission from reeonneetion zones. We diseuss their different eharaeter- 
isties and present a simple argument, whieh shows that in eontrast to “non-boosted” reeonneetion 
models (i.e. the ease where the partieles’ veloeities are not strongly an-isotropie in the jet frame), 
anisotropie models ean be eonsistent with the observed variability in GRB light-eurves. Seetionj^ 
is devoted to a detailed deseription of the model that is eonsidered in this work. Next, in § [^analytie 
expressions are derived for the resulting light-eurves, and their dependenee on the model parame¬ 
ters is diseussed. Seetion § [^features an extensive eomparison of our results to GRB observations. 
Finally, our main results are summarized and diseussed in § [^ 
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2. Motivation 


2.1. The Basic Setup for Relativistic Magnetic Reconnection 


In GRBs and other relativistic outflow sources, magnetic reconnection is considered to be 
potentially most relevant if the outflow is initially Poynting flux dominated. In this case magnetic 
reconnection is expected to be most efficient and prominent if a high a is maintained in the outflow 
out to the reconnection site, so that the plasma that is inflowing into the reconnection region (or 
layer) is highly magnetized, with a S> 1. Such a situation is referred to as relativistic reconnec¬ 
tion, as the total energy largely exceeds the rest energy of the particles and the outflow from the 
reconnection region is expected to be relativistic. For such relativistic magnetic reconnection, two 
main different types of reconnection models have been discussed in the literature. 


The first type involves quasi-steady-state reconnection configurations with a continuous in¬ 
flow of magnetized plasma into the reconnection layer, and a continuous outflow out of it. Analytic 
models of relativistic magnetic reconnection (in which the inflowing plasma is highly magnetized 
with cr 3> 1, e.g., Lyubarsk^ 20051 assume a single velocity for the plasma out-flowing from 
the reconnection region, which they find to be relativistic (corresponding to a Lorentz factor of 
r' ~ 1), and do not distinguish between individual particles of different energy. Particle 


In Cell (PIC) simulations (e.g., [Cerutti et al. ||2012[|2OT^|Sironi & Spitkovsky ||2014HGuo et al. 
2014[ Melzani et al. 2014[ Kagan et al. 20151, on the other hand, find that electrons that are ac¬ 
celerated to higher energies also spend a longer time being accelerated in the electric field within 
the reconnection layer, near the x-points, and thus their velocities become more tightly collimated. 
This leads to a positive correlation between the random Lorentz factor of the electrons 7 ' and their 
bulk Lorentz factor relative to the jet frame, T'. Moreover, once the electrons are deflected out 
of the x-point they quickly enter magnetic islands where their velocities are isotropized over their 
Larmour gyration time. Therefore, only electrons that emit a good fraction of their radiation within 
less than their Larmour gyration time will produce highly beamed radiation in the jet’s bulk frame. 
This corresponds to emitted photon energies above the synchrotron burn-off limit (> 100 MeV in 
the jet’s bulk frame), which is also known as the “maximum synchrotron energy” ( de Jager||1996 ) 
and corresponds to an observed energy of hi'syn,ma.^ = -^syn.max ~ 7(1 -f ^)“^(r/100) GeV. Such 
high-energy electrons near the burn-off limit are highly ’’fast cooling”. We stress that this set-up is 
therefore likely less relevant for explaining properties of the < MeV prompt emission. Within this 
model, the low energy part of the spectrum is likely to be produced by non-boosted synchrotron 
emission in the jet’s frame. Since the same electrons are producing both anisotropic and non- 
boosted emissions, the spectrum is expected to transition smoothly from the former to the latter, as 
the observed frequency decreases. Various implications of this scenario (at higher frequencies) are 
considered in detail in 


4.5 


In this picture, in the vicinity of the x-points where > B'^ the electric field cannot vanish in 
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any rest frame, and the partieles are direetly aeeelerated in the eleetrie field. In the magnetie islands 
j ^2 ^ j ^2 partiele veloeities have been isotropized sueh that the eenter of momentum (CM) 

partiele veloeity is the same for partieles of all energies, so that in this plasma’s rest frame, 
the eharged partieles ean short-out the eleetrie field, eausing it to praetieally vanish (and thus 
approaehing the ideal-MHD limit within the magnetie islands). As the veloeity of the islands 
relative to the jet’s bulk rest frame is small (/^j^i 1 ) the eleetrie field in the jet’s bulk frame 

is also very small there {E'/B' ~ 1). This may be assumed to still hold to zeroth order 

also in the intermediate region between the x-points and the magnetie islands, where most of 
the beamed radiation is expeeted to be emitted. In this region one already has but 

the partiele veloeities are still highly anisotropie and the partiele CM veloeity is energy dependent, 
f^CM = /^cm( 7D’ that anyway the eleetrie field eould at most vanish in the CM frame of partieles 
of one partieular energy but not in that of partieles of all other energies. 


The seeond type of relativistie magnetie reeonneetion models invoke relativistie turbulenee 


(|Lyutikov & Blandford 

2003 

Kumar & Narayan |2009 

Lazar et al. 

2009; 

Inoue et al. 

2011 ; 

Zrake 20 14^ 

2015 

; Lazarian et al. 

20151 East et al. |2015t Zrake & East|2016). In the eontext of 


GRBs within this type of model, different regions of the jet undergo sporadie reeonneetion events 
that give rise to relativistie bulk motions of plasma - “blobs” or “plasmoids” - relative to the jet’s 
frame. In this pieture eaeh blob is maeroseopie and signifies a relativistie bulk motion relative to 
the jet’s frame, where the blob’s bulk veloeity is equal to eleetrons of any 7 '. 


Notiee that although many plasmoids are found to be generated also in reeent PIC simulations 


( |Sironi & Spitkovsk^|2014j [Sironi et al. ||2015| ), this strueture is eontinuously evolving in time, 
as the magnetie islands merge with eaeh other and grow larger and larger. This situation eould 
then lead to an evolution of the speetrum within the prompt phase that is not observed. However, 
if the emission from all of the magnetie islands eomposing eaeh reeonneetion layer produees a 
single spike in the light-eurve then it would prediet speetral evolution within a single spike, whieh 
is observed. Moreover, sinee the inerease in islands’ sizes during this evolution is roughly linear, 
this would also imply that the aeeeleration time-seale of the highest energy eleetrons is dominated 
by the growth time of the largest island, whieh might explain the late onset of the high-energy 


emission (eompared to the < MeV emission) in most bright Fermi/LAT GRBs (e.g., Abdo et al. 


2009a|a|bt [Aekermann et al.||2010[ ). However, in this pieture during a signifieant fraetion of the 
pulse’s dynamieal time, it is no longer possible to maintain face ^ fdyn as required in order to have 
eontinuous heating that balanees the radiative eooling in order to avoid exeess optieal and X-ray 
emissions. Furthermore, reeent studies ( Melzani et al.||2014 | find that the plasmoids generated in 
this setup are not neeessarily traveling relativistieally with respeet to the jet’s bulk frame. A blob¬ 
like emission deseribed by the k = 0 mode is, however, observed in numerieal studies involving 


relativistie turbulenee. In this ease, the turbulenee ean enhanee the reeonneetion rate (e.g., Lazarian 
etal.|2015t[Zrake|2015D . 
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Ideal MHD holds within each blob and the electric field within its volume vanishes in its own 
frame. The electrons maintain a roughly isotropic velocity distribution within the blob’s rest frame 
and their emission is approximately isotropic in this frame. Hence, the blob’s emission is beamed 
in the jet’s rest frame because of the blobs relativistic motion in this frame. Previous studies have 
focused on turbulence that is isotropic in the jet’s frame. However, given an initially striped wind 
configuration of the magnetic field, it is likely that this turbulence would be largely constrained to 
a thin layer parallel to the shell front. Such an anisotropic relativistic turbulence would result in a 
correspondingly anisotropic emission in the jet’s frame. 


In this scenario the electrons may be slow cooling. Moreover, second order Fermi acceler¬ 
ation within the turbulent region can continuously heat the electrons until they reach a balance 
between heating and cooling (e.g., Murase et al. 2012; Asano & Terasawa 2015| ). This naturally 
allows for the electrons to remain hot while radiating, which is necessary in order to explain the 
observed prompt spectra from GRBs in a magnetically dominated emitting region (Beniamini & 
Piran ||2014 ). In addition, boosted emission is necessary in order to explain the observed prompt 
variability as will be shown below. For these reasons, in this work we focus mainly on this second 
scenario, while the first scenario is addressed mainly in 


4.5 


Finally, we remark that although the A; = 0 and k = 1 models involve different modes of 
reconnection, they both rely on a similar initial physical set-up: a magnetically dominated jet, 
where each pulse in the light-curve is produced by some macroscopic reconnection event in a layer 
perpendicular to the jet’s propagation direction. While the fc = 0 is more likely for explaining the 
observed properties of the sub-MeV pulse, it is possible that in at least some GRBs, the k = 1 mode 
is responsible for producing the GeV emission. In this case, the two modes may be operating in 
the same outflow. This may occur, for instance, if the initially quasi-steady state eventually drives 


turbulent reconnection until the later becomes the dominant mode of reconnection (see |Lazarian et 
aH2015 ). The ratio between the GeV emission and the sub-MeV pulse would depend in this case 
on the fraction of the flow involved in each reconnection mode and on their respective durations. 
Whether or not this can reproduce the observed ratio of these spectral components could be a 
potential test for this scenario. 


2.2. Lightcurve Variability 


In a simple, non-boosted emission picture, the reconnection model cannot easily reproduce 
the observed variability. This is because in order to obtain large-scale ordered reconnection, the 
speed of the incoming plasma that approaches the reconnection region from both sides is expected 
to be v[^ = I3[^c < 0.1 — 0.25c ( Lyubarsky]||2005| Guo et al.|2015 ; Liu et al. '||2015 1. Therefore, 
if the radial width of the region feeding the reconnection layer is L in the lab frame and L' = TL 
in the jet’s bulk frame, then the reconnection lasts At' = L'/v'^^ in the jet’s bulk frame, and 
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At = r At' in the lab frame. The time over whieh the radiation that emitted from this reeonneetion 
event along the line of sight reaehes the observer, or the radial time of the resulting pulse, is 
ATr = (1 — (3) At ~ L/2v[^. The inerease in radius during the reeonneetion event is AR cAt, 
and sinee this is a lower limit to the emission radius, R > AR, the angular spreading time in the 
arrival of photons to the observer, whieh is a lower-limit to the pulse width, is at least 


AT« = 


R 


> 


AR 


2cr2 - 2cr2 2vi 


ATr 


( 1 ) 


The total pulse width aeeounts for both radial and angular spreading in the photon arrival times, 
and is therefore 

AT ~ ATj, + ATg ~ —— . (2) 

This is a faetor of ~ 10(/3'j^/0.1)~^ larger than the differenee in the ejeetion time of neighbouring 
reeonneeting layers or shells, Afgj ~ T/c, whieh would also be the differenee in the onset times 
of neighbouring pulses for emission from the same radius. While fluetuations are possible around 
this typieal differenee in pulse onset times, the mean separation does not vary unless the emis¬ 
sion radius systematieally ehanges throughout the GRB. There is, however, observational evidenee 
against this, sinee the gross properties of the GRB emission do not appear to systematieally ehange 
throughout the prompt emission. 


The time between pulses ean inerease if only a small portion of the jet material would have 
eontributed to the observed radiation (i.e. if the distanee between the reeonneetion layers is larger 
than L, so that Afej would be eorrespondingly larger, where the added outflow layer does not 
eontribute to the reeonneetion). However, assuming a roughly eonstant jet power this would imply 
that only a eorrespondingly small fraetion is dissipated and eontributes to the gamma-ray emission. 
This is beeause signifieant pressure gradients tend to be washed out, so that a smooth pressure 
profile may be expeeted in the propagation direetion. For a fixed loeal bulk Lorentz faetor the 
energy flux is proportional to the proper enthalpy density, w + B'^/An = w(l -f a), whieh in 
the highly magnetized regions {a S> 1) is dominated by the magnetie term, half of whieh arises 
from the magnetie pressure B‘^/8n. The thermal pressure, on the other hand, ean aeeount for 
at most one quarter of the partiele proper enthalpy density, w = pc^ + Cint + p, sinee even for a 
relativistieally hot plasma p = Cint/S ~ w/d. Therefore, a eomparable energy per unit length in the 
propagation direetion is expeeted also for low magnetization regions (a < 1) or highly magnetized 
regions that do not border an oppositely oriented magnetie field region with whieh they ean easily 
reeonneet (e.g. if the field orientation ehanges randomly due to fluetuations in an aeeretion dise 
around a blaek hole, rather than very orderly as in a strong MHD wind launehed by a newly born 
milliseeond magnetar). This eould be very ehallenging for the total energy budget of the GRB, as 
the total jet energy would be even larger than eurrent in estimates, whieh are already eonsidered 
highly demanding for viable GRB progenitors ( Panaiteseu and Kumar||20d2j [Granot et al. 2006 
Fan et al.|2006i|Be^niamini et al. |2015] ). 
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All this implies that the “simple” reeonneetion model prediets pulses that are significantly 
broader than the intervals between them, contrary to what is seen in observations (see also [Lazar 


et al. 2009), namely that these times are roughly equal (jNakar & Piran 2002). Therefore, it is 


interesting to explore a slightly more complicated yet physically motivated reconnection model, 
which may be able to both account for the continuous heating, as well as reproduce the observed 
variability. In this model the plasma in the thin reconnection region is ejected away from the X- 
point along the reconnection layer at a relativistic speed in the jet’s bulk frame, corresponding to a 
Lorentz factor T'. The emission is assumed to be isotropic in a frame moving at T' relative to the 
jet’s bulk frame, and is therefore anisotropic in the jet’s bulk frame for T' S> 1. 


Anisotropic emission in the jet’s bulk frame only reduces ATq by a factor of ~ T', while it 
does not change AT^. Therefore, given that AT ~ AT^ + ATe, for a smooth continuous emission 
as a function of radius (and T' > a few) AT ~ AT^ ~ R/2cT‘^ Afej/2/3i'n = 5(/)jn/0.1)“^Afej- 
However, since ATq ~ AT^/T' is significantly smaller than AT ~ AT^, and the emission at any 
given radius originates from a small part (a fraction ~ 1/(T')^) of the visible region (corresponding 
to angles 6 < 1/T from the line of sight), variations in the emission intensity with radius would 
be reflected in the light curve down to the shorter angular time-scale AT^/T'. Moreover, such 
variations are arguably quite likely in efficient (high a) magnetic reconnection (i.e. reconnection 
leading to dissipation of a significant fraction of the total energy), which is rarely very smooth and 
continuous. Therefore, it is likely to typically observe variability on the shorter angular time-scale 
ATr/V. 


One implication of this model is that it naturally results in larger emitting radii than those 
usually implied from the variability time-scale. This is because for non-boosted emission in the 
jet’s bulk frame the observed variability time-scale is given by: 


A ry-i _ R'lSO 

iso — 


2cr2 ’ 


(3) 


whereas in the current model it is shorter by a factor of ~ T', AT ~ ATiso/T'. This would 
imply a larger emission radius in order to reproduce the same observed variability time-scales. For 
r 300, T' 10 and AT = Isec we obtain i? ~ 5 x lO^^cm. This radius is ten times larger than 
would be inferred for a non-boosted emission in the jet’s bulk frame, and is still a few times smaller 
than typical estimates for the deceleration radius (and therefore consistent with the interpretation 
that the prompt signal originates at R). 
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3. The model 


3.1. The physical model 


We associate each pulse in the light-curve with emission from a different thin “shell”. In the 
context of a striped wind each “shell” is associated with the thin reconnection layer between regions 
(wider spherical shells) of oppositely oriented tangential (primarily toroidal) magnetic field lines. 
In this section we discuss the emission from a single shell, and in § 4.6 we show that adding up the 
emission from many shells produce light curve variability similar to that seen in observations. 




Fig. 1.— Left: Schematic geometry of the model at hand. The emitting matter is moving relativis- 
tically out of the reconnection x-points at a Lorentz factor F' in the jet’s bulk frame (i.e. relative 
to the bulk of the outflow). Right: The two reference frames used in this work: the lab frame 
with the origin at the central source (in either Cartesian or spherical coordinates), and the jet’s bulk 
frame (denoted by primes), travelling at a Lorentz factor F in the f = z' direction and centred at 
the emission point. 


Consider an ultra-relativistic (moving with F S> 1 in the rest frame of the central source) thin 
(of radial width ^ Rj (F^F')) “spherical” (occupying an angle Q around the line of sight such that 
> a few) shell which starts emitting at some initial radius i?o (possibly, but not necessarily due 
to reconnection) and turns off at a final radius Rj > i?o- Let the emitters be moving with a Lorentz 
factor F' in the local bulk frame of the flow, in two opposite directions that are perpendicular to the 
radial direction. There are three relevant frames in this paper: (i) the central source frame, (ii) the 
jet’s local bulk frame, which moves with a Lorentz factor F in the radial direction to the central 
source, and (iii) the emitters’ frame which moves with a Lorentz factor F' compared to the jet’s bulk 
frame, in which the radiation is assumed to be emitted isotropically. Quantities in frame (i) are un¬ 
primed, those in frame (ii) are primed, while quantities in frame (iii) are double primed. With these 
definitions, the emitters’ speed is given by: jS' = ±/3'x' {z' = f is the radial direction in the jet’s 
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bulk frame), x' is the projeetion of x (the direetion perpendieular to the direetion to the observer 
n = z) on the plane perpendieular to the radial direetion. Notiee that: n' = f cos 9' — 9 sin 9'. In 
addition: 


{9 ■ x)9 + {(() ■ x)<:f) 0 cos 0 cos 0 — 0 sin 0 


^(0 . £)2 . ^^2 \/sin^ 0 + cos^ 9 cos2 0 

9 cos 9 cos 4) — 4) sin (j) - - 

9 cos (p — (p sin 4 >, 


(4) 


\/l — sin^ 9 cos^ 4> 

where /5' = ±x' and we have approximated cos^ 9 = 1 — sin^ 0 1. Moreover, 

cos 9 cos 4> 


cos 4> = x' ■ 9 = 


\/l — sin^ 9 cos^ 4> 


COS0 


(5) 


Combining this with the expression for n' we get that: n' • /3' = ± sin 9' cos 4>' ■ For a schematic 
figure of the geometry see Fig. We aim to find the resulting light-curve in the observer’s frame. 
Our modelling follows that of [Genet &: Granof] (|2009|) (which relies on the formalism developed 


m 


Granot | ([20051) and|Granot et al.|([20081)), and converges to the results in this paper in the limit 


of r' —)■ 1 (and (5' —)■ 0), corresponding to the case of non-boosted emission. Notice that once the 
motion of the local center of momentum velocity of the emitting particles becomes non-relativistic, 
their radiation is no longer beamed and even if they emit in a very ordered magnetic field, their 


emission pattern will be broadly similar to isotropic emission in the jet’s bulk frame (Granot et al. 


2002 [ [Granot fc Konigl 2003 [ Beloborodov et al. 2011 1 , and the situation reduces back to the 
simple “one zone” case. 

The jet may be accelerating, coasting or decelerating as it is emitting (see § [^for more details). 
To account for this, the Lorentz factor of the bulk is assumed to evolve as a power-law with radius: 

= Tl{R/R q)~^ where Fq = r(i?o)- The case of a coasting jet (with no acceleration or 
deceleration) can be easily obtained by setting m = 0 in the following equations. The luminosity 
is assumed to be isotropic in the emitter’s frame (double primed frame) and depends only on the 
radius: 

In the last expression, f{x) is a dimensionless function of the radius that is normalized such that 
f{x)dx = 1, z/" is the frequency where peaks, and S{x) is the spectrum (which is 

assumed to be a Band function). 


Six) = . 


(7) 


where Xh = (as — /5 b)/(2-|-q;b), while as and /3 b are the high and low energy slopes of the 
spectrum. For as > —2 > /3 b the Band function has a peak in the uFy spectrum, at x = 1, 
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and therefore sinee S{x) is normalized sueh that S{x) = xS{x) = 1 at x = 1, it will not affect 
normalization of uFi, at its peak. The two functional forms used in the Band function are matched 

at = XbJy'p. 

The flux is then given by: 


FAT) 


1 

AttD^ 


dL^ 


1 

47rT)2 


V^dL'^, 


1 

AttD^ 


'ry/3-k 

V^^—dL'A 


( 8 ) 


where for convenience we denote the effective distance to the source in terms of the luminosity 
distance and red-shift hy D = di/y/l + z. In addition, V is the Doppler factor from the jet’s bulk 
frame (the rest frame of the bulk) to the observer frame (the rest frame of the central source), V' is 
the Doppler factor from the emitter’s frame to the jet’s bulk frame and dL”,, = L'l„ {R)dnd(l)/ (dvr). 
In the last expression, /i = cos 9 where 9 is the angle between the emitting location and line of 
sight and 0 is the azimuthal angle (these are the angles in a spherical coordinate system whose 
origin is at the central source while its z-axis points toward the observer). This implies that the 
observed frequency (u) is related to the frequencies in the jet’s bulk frame (z/') and in the emitter’s 
frame (u") by (1 + z)u = V u' = W'v”. The parameter k can obtain only two possible values: 
/c = 1 in case each emitter produces continuous emission (a steady state in the jet’s bulk frame, 
e.g. as in 14.5), or fc = 0 in case the emitters produce a single “blob” (as is likely for the turbulent 
scenario discussed in Q. The observed time T is set such that T = 0 corresponds to the arrival 
time of a photon emitted at the central source (i.e. in the origin) at the ejection time of the shell. 
The first photons to reach the observer from the emission at radii Rq < R < Rf are from along 
the line of sight {9 = 0) from R = Rq, and arrive (accounting for cosmological red-shift) at 


(1 -f z)Rq 
2{m + IjcTp 


(9) 


We define the equal arrival time surface (EATS) as the locus of points from which photons emitted 
at different radii R and angles 9 relative to the line of sight, reach the observer simultaneously at 
an observed time T. The maximal radius along the EATS, Rl{T) (which is obtained along the line 
of sight), and the Eorentz factor at that radius, Tl{T), can be conveniently expressed in terms of 
this time. 


^L_fRL\ ^ _fT\ 

l^o~\To) ’ ~\To) 


( 10 ) 


3.2. The observed flux 

In order to calculate the flux density F^, that reaches the observer at time T we integrate the 
luminosity L"„ over the EATS. This introduces a relation between R and 9,T/{l + z) = t — Rjji/c, 
for photons arriving at a given time. In addition, since the emission in our model is anisotropic. 
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photons emitted at different 0 will have different Doppler faetors from the emitter’s frame to the 
bulk frame and we eannot use symmetry arguments to integrate over 0. As we assume F 3> 1, 
beeause of relativistie beaming only angles of 6^ < a few eontribute signifieantly to the observed 
flux, so one ean safely assume 6^ ^ 1. It is eonvenient to use the dimensionless normalized 
radius y = R/Rl and ehange variables of integration for ealeulating the flux from y io y ox 
dy —)■ \dy/dy\dy, using the relation 


dy 


dy 


y + my"^ ^ 


T \ ""+1 1 (m + y ^ ^)y'^ ^ 


2(m + l)r| \To J Fq 2(m + 1) 


( 11 ) 


The Doppler faetors and ean be expressed as 


V = 


V = 


V{l-Py) 


= Ti 


2{m + l)y-^R f T \ 2(m + l)y-^R 


m + 

sin 6' = 2 


TqJ " m + y 


—m—l ’ 


F'(l —/9'sin6*'cos0') ’ m + y~'^~^ 

One ean also define z/q = 2FoFVq /{I + z) and use it in order to express 

m 

v” V / ■?/ \ 2 


( 12 ) 

(13) 




(■ 


m + 1 


y 

l/mir 


F'^(l—/9'sin6*'cos(/)) , 


(14) 


where the limits of integration over y are from y m\n = min(l, Rq/Rl) to z/max = min(l, Rj/Rl). 
A general expression for the flux is then given by: 


2FoF'L^„ / X \ 2(m + l) rj/max 

MT)=—r^T^{^] I dy 


AttD^ \ Tr 


' J/mii 


m + 1 


/—m—l 


y 


_2_nL 


/ 


y 


2 "' \ rn + l 

fn 


X 


27rF'4 


m+y~ 

p27r 

I (i(/) (1 —/?'sin6*'cos 0)^“^ S'[x(0, ?/)] . 


(15) 


4. Results for single-pulse emission 


We turn to a detailed diseussion of the results of our model. In §4.1[ 4.2 we present simple an- 
alytie expressions for the general flux provided in Eq. [^for speeifie speetral and radial emission 
profiles. The light-eurves of single pulses under various assumptions are presented in (4.3 Their 
dependenee on the radial width of the jet eontributing to the emission and on the emitters’ Lorentz 
faetor, F', as well as estimates for their degree of asymmetry are taken up in (4.4 In (4.5 we 


deviate from the main model eonsidered in the paper, and, motivated by PIC simulations, turn to 
a diseussion of the predieted emission from a eontinuous injeetion of partieles to the reeonneetion 
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layer resulting in a eorrelation between the random Lorentz faetor of the eleetrons and the bulk 
Lorentz faetor of the emitters. The results derived under these assumptions are likely to be less 
relevant for explaining the properties of the sub-MeV pulse, but ean still have interesting impliea- 
tions for the emission at very high frequeneies. Finally in §4.6[ we present the results for the entire 
light-eurve using different distributions of the intrinsie parameters of the model. 


4.1. analytic estimates for power law spectra and power law function in radius 


For a power law speetrum, S'pl(x) = one ean obtain a somewhat simpler 

analytie formula for the observed flux. In this ease sinee f{R/Ro) = 0 for R < Rq (see below), 
one ean effeetively plug j/mm = in Eq. ([T4|) to obtain: 


v” V / T \ 2(™+l) ( m + y^ rr, . . 

^ = (l-^^smeco.^), 


(16) 


Substituting S{x) —)■ S'pl(x) = x " and: 


f(R/R^) = J Ro< R< Rf , 

0 otherwise, 


(17) 


in Eq. ( fTSj ) yields 


MT) = 


2TorLl 

AttD^ 




2a — m{l-\-ct) 

R \ 2(m + l) 


(a + 1) 


^ 2/max 


dy 


m+1 


m+y 


—m—1 


2+a 


y 


a-l-fil+a) 


X 


2 ?/V _ 1 Jymin 

Ro J 

1 

(1-/3'sin 0' cos 0)^-3-“ 


27rr'4+2« 


( 18 ) 


It is easy to see that Eq. (18) reduees to Eq. 13 of Granot et al. 


for r' = 1 (and without 

the attenuation term in that paper, due to 7-7 opaeity). At the other extreme, F' ::§> 1, Eq. (18 1 ean 
be approximately estimated (see Appendix]^ for a derivation) as: 


F.{T) 


-i-> -i->/ 1 T n 2a —777.(1 + 0 :) 777(0 —1)— 2(a+l) 

Tor / z/ \ / T \ 2 {m + l) (q _|_ X) (/Ti + 2) 2 (m+l) F (2 + Of 

2'=7rE)2 r(3 + a 


XI 

k) • 

(19) 


We see that for a blob like emission {k = 0) we obtain FuiT) cx F' as expeeted from beaming. For 
eontinuous emission {k = 1), the flux appears to be suppressed by a faetor F'“^. This is beeause 
in eontinuous emission, the emission at eaeh loeation in the jet’s bulk frame remains eonstant in 
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time as the emitting plasma flows through it, while for a single blob of plasma the loeation of the 
emission (i.e. of the blob) moves in the jet’s bulk frame and ehases after the photons it emits, and 
for the viewing angle from whieh it appears most luminous it trails them just by a small amount, 
so that the same emitted photons (over the same time interval in the jet’s bulk frame) arrive to a 
stationary observer (in the jet’s bulk frame, and also to us) over a shorter time interval (by a faetor 
of,KJdt' = (T'V)-' = 1 — /S' sin 9' cos 0') eompared to a steady, eontinuous flow, resulting in a 
flux larger by the inverse of this faetor, T'V. Alternatively, a eontinuous flow ean be approximated 
as eomposed of a large number of blobs, in whieh ease only a fraetion (r''D')“^,of the N blobs 
that are instantaneously within a given segment of the flow eontribute to the observed emission 
(or photon front) from that segment, resulting in a flux smaller by this faetor eompared to a single 
blob with the eombined luminosity of the N blobs. For the observed flux in our ease, the viewing 
angle relative to the direetion of the emitting fluid in the jet’s bulk frame ehanges with time and 
between the different parts of the souree, so that altogether the total flux refleet a weighted mean 
of whieh is ~ F'”^, for the suppression faetor of a eontinuous flow eompared to a single 

blob. Stated differently, the total energy emitted in the two oases, soales as E oc sinoe 

in order to emit over the same observed duration, a blob needs to emit in its own rest frame over a 
time longer by a faetor of ~ F'^ with the same luminosity. This implies that in order to keep the 
total emitted energy eonstant in the two soenarios, we must ohange the luminosity in the emitters’ 
frame aeeording to: L'l,, oc F'^^. Plugging this baek to Eq. (19) we obtain the expeoted sealing 
with F', independent of k. 


4.2. analytic estimates for power law spectra and single emission radius 


Another simple test-ease is that of a power-law speetra (as in (4.1) but with a delta-funetion 
emission in radius. 


/ 



Ro6{R-Ro) 




ymin^iy 2/min) ; 


( 20 ) 


where = (T/Tq) V(m+i)_ in this ease, one ean analytieally perform the integral over y and 
obtain: 


KiT) 



It is eonvenient to use here the initial radial time Tq = Tr{Ro) given by Eq. (|^, whieh is 
the time when the first photons reaeh the observer (for our definition of T = 0 as the arrival 
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time of photons from the souree when the shell is ejeeted), and the eorresponding angular time, 
Tg = TelRo) = Ro/{2crl) = {m + l)Tr{Ro) = (m + l)To that determines the width of the pulse. 
We also define Tg = Tq — Tg = —mTo. We may now write the expression for the flux as: 


Fu{T) = 


/lyy- / t-t, 

[uoj V Tg 


-2-a 


■p/—4—2a r2TT 


-2(5' 


V(T-To)T, 

T-Tg 



k—S—Oi 


( 22 ) 


For (5' = 0 the integrand beeomes 1 and we obtain a power law deeay with a temporal index 
2 + a as is well known for high latitude emission. At the other extreme, for /)'—)■ 1 the integrand 
beeomes maximal when sin 9' = 2 '^/(rTTTTjT/(m + 1 + T) 1 where T = {T — Tq)/Tq, whieh 
oeeurs for T m + 1, and cos0 ~ 1, whieh oeeurs for 0 0. Therefore, the peak in the light 

eurve oeeurs at small values of f = [T/(m + 1)] — 1, and for estimating the fluenee (i.e. evaluating 
the light-eurve around its peak) we may assume that f <C 1, where in this limit sin^' 1 — ?/8, 
and (T — Tg)/Tg = (m + T /Tq) /(m + 1) ^2. In terms of 0, the width of the region where the 
dominant eontribution eomes from is roughly ~ l/V -C 1 around 0 = 0 or l^jF' < a few, so that 
we ean safely assume |0| ^ 1, cos0 ~ 1 — 0^/2, and take the limits of integration to be from 
—oo to cx), where 1 — /)'sin6*'cos0 [1 + (r'f/2)^ + (r'0)^]/2r'^ = (! + ? + (jF)/2V'‘^ (where 

i = r'f/2 and 0 = r'0). 


R. 


[T \ FTqT'L'' / 1/ \ {T-Tg 

TqT'‘^~‘^^L"■■ / \ roo 


2\k-3-a 


^ /.y° r(2.5 + a-fc) 

2^-i-“7r3/2_D2 yi/Qy r(3 + a —/c) 




(23) 


We have dT = {2Tg/T')dt, and near the peak F^{t) ^ K{1 + ^ ^ that the fluenee is given 

by 




dTF,{T)^‘^-^ rdi{l+P)’‘-l- 

^ J —OO 


/^\- r(2 + a-fc) 

2 k - 2 - a ^£)2 p (3 _|_ q , _ / j ) 

(24) 


The dependenee on Tq and T' in Eq. (24) is the same as in Eq. (19), while in Eq. (23) there is an 
extra power of T'. This arises sinee the width of the peak in the light-eurve for the emission from 
a single radius is ~ Tg/V, i.e. ~ l/V 1 times the angular time from that radius, due to the 
beaming of the emission in the jet’s bulk frame. The time dependenee of the emission near the 
peak, in Eq. ([23]) is dominated by the beaming in the jet’s bulk frame, while at late times T ^ Tg 


15 





















it becomes dominated by the usual high-latitude emission, as the beaming in the jet’s bulk frame 
approaches a constant (i.e. the n' — f and the integral in Eq. ( |2^ approaches a constant, leaving 
only the explicit time dependence in front of the integral). While for AR = Rf — Rq ~ i?o the 
pulse width is ~ Tq, the normalization of the emission from a single radius is such that the fluence 
in both cases is similar. 




Fig. 2. — Lightcurves of a single pulse with a power-law between Ro and Rf = 2Rq {left panel) 
and log-normal {right panel) dependence of the luminosity on radius. All cases are plotted for a 
typical Band function spectral emissivity and E' = 3, m = 0, /c = 0. The flux is given in units 
of the “typical” peak flux Fq (giving in Eq. p6|), and shown as a function of the normalized time 


T or T (given in Eq. [271). T = 0 corresponds to the ejection time, and T = 0 corresponds 
to the time when the first photons reach the observer. We show results for power-law and log¬ 
normal luminosities as a function of radii for three different values of the power-law index a and 
the standard deviation ainR, respectively. 


4.3. Band spectrum emission for different dependencies of the luminosity on radius 


Using Eq. ( 151 we calculate the resulting light-curve for a Band spectrum emission function 
with a low-energy photon index as = —1 and a high-energy photon index /3s = —2.3, typical 
values obtained from fitting of GRB data. For the radial dependence of the emissivity we use some 
specific functions for /{R/Rf}. For clarity, in all cases we take F' = 3, m = 0, A; = 0. We plot two 
typical cases in Fig.[^ a power law in radius (Eq. [Tt]) with Rf = 2Rq and a Gaussian in ln(i?): 


f{R/Ro) = 


Rq! R 


exp 


[ln(i?//(o,eff)]- 

2crinR^ 


(25) 
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where -Ro.efr = exp((Tinj?^)-Ro is a slightly different effeetive radius ehosen sueh that /{R/Rq) 
would peak at i? = Rq. Note that ai^R can be related to the typical width in R through: AR 
sinh((JinK)-Ro (which for ainR 1 yields AR ^ cxinij-Ro)- We plot the results at the frequency Up. 
The flux is given in units of the “typical” peak flux, 


2ror'i-2^L" 
— _ ^ 


(26) 


and the time in units of the typical time Tq, either starting at the observed ejection time Tej (denoted 
as T) or starting at Tej + Tq, the time when the first photons reach the observer (this is defined as 

f), 

To 


T = 


T 


T ■ 


T = 


T 


T • 


(27) 


To To 

For a power-law emissivity, f{R/Ro) (x the pulse shape becomes less symmetric and peaks 
at earlier times as a decreases. For a Gaussian f{R/Ro), on the other hand, the peak occurs later 
for smaller typical widths (AR) but becomes more symmetric. The power-law and Gaussian cases 
approach the case of a delta function in radius and reproduces the result from §4.2| in the limits 
a ——cx) and AR —0, respectively. We note also that in the Gaussian case the typical width of 
the pulse is ~ max(To/r', ToAR/Ro) in accordance with Eqs. ( |2T[ ), ( [T^ . 




Fig. 3.— Left: single pulse light-curves for different assumptions on the emission (continuous 
and steady in the jet’s bulk frame vs. single blob vs. non-boosted emission). All cases are plotted 
with a typical Band function spectral emissivity, a constant dependence of luminosity on the radius 
with Rf = l.li?o as well as m = 0. The flux is given in units of the “typical” peak flux Fq as a 
function of the normalized time T. Right: light-curves with fc = 0, F' = 3, a constant dependence 
of luminosity on the radius with Rf = 1.1 To and different values of m. 
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4.4. Exploring the single-pulse parameter space 


Here we explore the dependenee of the single pulse light-eurve on the various model parame¬ 
ters. We define {T^y) as the first time when the flux first (last) reaehes Y% of its peak value. 
The rise (decay) time, is then defined by: Tyse.x = -Tr,x_ (Tdecay.x = Td,x_ -Td,x+) where 
X± = (100 ± X)l2. In addition, we define the “central” time, Tcy = — ^r,x+ (see FigQ. 

Notice that if there is a significant “plateau” stage then it is possible that Tj-ise, Tdecay ^ T^. If the 
rise and decay of a pulse are monotonic, then Tioo-x. = Tnse,x+Tcy+Tdecs.y,x = Tdy_ —Try_ is 
the time spent at a flux above X_% of the peak flux (analogous to the definition of GRB durations 
but with flux instead of fluence). 


^peak 



decay,50 


Fig. 4.— Schematic figure of the rise, decay and “central” times defined in the text for X = 50. 
Trise,x (Fdecay.x) IS the time it takes the flux to rise (decay) from X_% (X+%) to X+% (X_%) of 
the peak value (where X± = (100 ± X)/2). Tcy is the time the flux stays above X+% of its peak 
value. 


First, let us consider the more familiar case of non-boosted emission in the jet’s bulk frame. 
If XR < Rq then the flux peaks at approximately: T = Tgj + Tf, where Tf = {Rf /Rq)'^~^^To, or 
Tf = (Rf/Ro)'^'^^ — 1 = [1 + {AR/Rq)]'^~^^ — 1. The rise time is then given by: 

m+1 


T 

ri 


Tn 


(+ 

m+1 

r AR] 

- 1 = 

1 + c. 

\Ro/ 

1 

Rq _ 


- 1 


ar<s:Ro 


» (m + 1) 


AR 

Rn 


(28) 


and the decay time Tdecay ~ Te{Rf) is given by: 

AR 


ToJ ’ 

(29) 

leading to a total pulse width AT = + Tdecay of: 
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which for m = 0 leads to AT/Tq ~ Rj/Rq. For Ai? 3> Rq both the rise time and the decay 
time may be altered due to the emissivity funetion (see Table for details), whieh (for the ease 
of power law luminosity as a funetion of radius and frequeney) dominates the emission so long as 
0.5 < T < Tf. The relevant parameter in this ease is the power law index q in the expression for 


FyiT) cx T*?. Using Eq. 19 we define: 


0 . = 


2a — m{l + a) 
2(m + 1) 


(31) 


Speeifieally for q = 0 the flux at this interval reaehes an approximate plateau. This is an illustration 
of a situation in whieh 3> Rise- We define a quantitative measure of pulse asymmetry as A = 
Trise/Tdecay- Pulses are asymmetrie if A 1 or A ^ 1. In faet, even for A = 1, the shape of the 
light-eurve could still be different between the rise and deeay parts. For the non-boosted emission 
we find that unless the emissivity is inereasing with radius, then A 1 (estimates for individual 
eases are listed in Table implying shorter rise than deeay time and a large degree of asymmetry 
regardless of the relation between AR and Rq. In faet it is hard to produee very symmetrie pulses 
(A 1) for the case of such non-boosted emission, since the mechanisms controlling the rise and 
deeay times are different. We now eompare these results to the anisotropie ease. 



Fig. 5.— The emission region - the blue contours show the area within the emitting spherical shell 
that is within the beaming cone of the emitting material in the jet’s bulk frame, i.e. where n' is 
within an angle of l/T' from the direetion of its veloeity R = ±x', for T' = 3,4, 5, 7,10, 20. The 
big blaek eirele indieates an angle of 6^ = l/T around the line of sight to the eentral souree, while 
the red arrow shows the x'-axis. 


In Fig.j^we plot the dependenee of the pulse shape on the nature of the emission: eontinuous, 
single blob and non-boosted emission, for a eonstant eo-moving luminosity with AR/Rq = 0.1. 
For AR < Rq the flux peaks earlier for non-boosted emission than for the eases with anisotropy. 
The light-eurve is sharper and less symmetrie for the non-boosted emission so long as AR < 
Rq/V . This arises since for AR < Rq/V in the non-boosted ease the pulse rise time is determined 
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by AR (Trise ~ AR/Rq) while its decay time is Tdecay ~ Tq, whereas in the anisotropic case 
the pulse is rather symmetric (A 1), since Trise ~ ^decay ~ l/h' (and its width reflects that 
for the emission from a single radius). To easily understand the latter expression, consider that 
the emission from each radius is dominated by two regions of angular size ~ l/(rT) around 
the two points of the intersection of the x'-axis (where (3' = ±/3'x') and the circle 9 = l/T, or 
^ = (r^)^ = 1 (see Fig. 1^ and §& This dominant contribution arrives over a time To/r, 
which is thus the “smearing time-scale” compared to the very large T' limit. In the right panel of 
Fig. [^we plot also the dependence of the light-curve on m for the anisotropic case. As shown in 
^4.2| (and can be seen for the non-boosted case from Eq. 28 1 , in terms of the normalized time T 
the peak occurs at T ^ m + 1, which is when the beaming cone points toward the observer, after 
accounting for the aberration of light for the jet’s frame to the observer’s frame. This behaviour, of 
the peak occurring at a later time for a larger m, can indeed be seen in the right panel of Fig. 


In Fig. 1^ we plot the dependence of the pulse shape on the values of F'. Once F' 1 
and l/F' AR/Rq, the pulse’s shape does not depend strongly on F', and is set instead by 
/{R/Rq). This is because in this case, the typical times related to the emissivity function are 
larger than the pulse width from a single radius, AT ~ Tq/T'. The top-right panel of Fig. shows 
light-curves for AR/Rq = 1 and a log-normal radial dependence of the co-moving emissivity. It 
demonstrates that for Ai? ~ i?o with a smooth turn-on and turn-off of the emissivity, as might 
typically be expected, the pulse-shape for anisotropic emission is not very different from that 
of non-boosted emission, showing a similarly noticeable degree of asymmetry (A 50 ~ 0.107, 
Ago ~ 0.096 for the non-boosted case as opposed to A 50 ~ 0.118, Ago ~ 0.173 for the anisotropic 
case). The asymmetry becomes larger for AR S> Rq (and remains independent of F'). This 
is because the rise time is dominated by the radial time-scale at the point where the emissivity 
function is within la of the peak, i.e. Ri ~ R^/AR T^se ~ {Ro/^RJ^^, whereas the decay 
time is dominated by the radial time-scale where the emissivity function is la above the peak: 
Rf ^ AR Tdecay ~ (AR/Rq)"^^. TMs Icads to A ~ (Rq/AR^^'^^^ (in fact, there is a 
numerical pre-factor here of order 0.1 arising from the smooth shape of the log-normal function’s 
decline, see table [^. Slightly asymmetric pulses are in broad agreement with the typical observed 
pulse shape. On the other hand, symmetric pulses are also sometimes observed, and are hard 
to produce with non-boosted emission, but rather easy to produce for anisotropic emission. A 
summary of the typical times for the various cases is given in table The expected asymmetry 
measures in our model can be seen explicitly in Fig. |7] where we plot the value of A for different 
values of F', AR/Rq and for different assumptions on the emissivity as a function of radius. 
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Fig. 6.— Light-curves of a pulse assuming different values of the emitter’s Lorentz faetor F' 
(log 2 (r') = 0, 1, 2, 3, 4, 5) with m = 0 and fc = 0. The flux is given in units of Fq as a 
funetion of the normalized time T. Top panels: results for a log-normal dependenee of luminosity 
on radius with different typieal widths (ainR = 0.01, 0.1, 1 from left to right). Bottom panels: 
results for a eonstant dependenee of luminosity on radius with different typieal widths (AR/Rq = 
{Rf — Ro)/Ro = 0.01, 0.1, 1 from left to right). 
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K) 


K) 


case 

width 

r' 

f(R) 

Tri=e.X 

^decay,X 

A 

fa,X 

description 

1 

AR-^Ro 

r' = 1 

any 


TeiRf) 


Trise < ^decay 

asymmetric, no plateau 

2 

AR^Ro 

AR ^1^1 

any 

1 

P’ 

1 

P 

1 

Tc < ^rise + ^’decay 

symmetric, no plateau 

3 

AR^Ro 

1 < Afl < 1 

^ Ro 

PL; any q 

1 

P 

1 

P 

1 

~ ^ ^ ~ ^decay 

symmetric, short plateau 

4 

AR^Ro 

1 < Afi < 1 

log-normal 



1 

Rc < Trise + ^decay 

symmetric, no plateau 

5 

Afl>Ho 

r' = 1 

PL; g < 0 t 

Te(Ro) - 1 

Tf 



asymmetric, no plateau 

6 

Afl>Ho 

r' = 1 

PL; q = 0 

Te(Ro) - 1 

Te(Rf) 

¥ 

‘ 

As 

o 

T™«Ta«Tdecay«(^)™’^' 

asymmetric, with plateau 

7 

Afl>Ho 

r' = 1 

PL; q > 0 t 


fe(Rf) 

1 

771+1 

Tc < Trise + ^decay 

asymmetric, no plateau 

8 

Afl>Ho 

r' > 1 

PL; q < 0 t 

1 

P 

Tf 


Tc <C TVise + ^decay 

asymmetric, no plateau 

9 

Afl>Ho 

r' > 1 

PL; q = 0 

1 

P 

TeiRf) 

r' 

m+1 ^ ^ 

Tc W ( ^ » ^deoay » ^decay 

asymmetric, with plateau 

10 

Afl>Ho 

r' > 1 

PL; q > 0 t 


'i'eiRf) 

r' 

r' 

771+1 

Tc ^ T^ise » ^decay 

asymmetric, no plateau 

11 

Afl>Ho 

any F' 

log-normal 


/ \ m-f-l 

•l0T^»10(^j 

/ P X 2(m+l) 

o-i(5ft) «i 

T™<CTc<Tdecay 

asymmetric, no plateau 


^assuming that (100 — X)/2 < f(Tf)/f{TQ) < (100 + X)/2. In case/(T_^)//(To) > (100 +-X’)/2 the situation reduces back to that obtained at g = 0, while if/(T^; )//(f’o) < (100 — X)/2 then the decline time stops increasing 
with AR and instead remains fixed on the value ~ obtained for f{Tf )//(T'o) = (100 — X) j2 with A ~ (P^Ty* oc 1/V'. 

^assuming that/(T_^)//(To) > . In case the increase with radius is shallower, the situation reduces back to that obtained at q = 0. 

* as it turns out, due to the smooth shape of the log-normal function, the decay time is approximately 10 times larger than Ty. 


Table 1: Properties of observed light-eurves for single pulses with different assumptions on the Lorentz faetor of the emitters (T') and q the 
power law index of T in the expression for Fy{T). 




Aso, log-normal 


A 50 , PL (a— 1 ) 



r' r' 


A50, PL (a= 0 ) 



<1 


A50, PL (a=l) 



r' 


Fig. 7.— The asymmetry ratio: A 50 = Trise, 5 o/^decay, 5 o as a funetion of F' and AR/Rq. All plots 
are for A: = 0,m = 0 and a power law spectrum with a = 1. Top Left: results for log-normal 
emissivity as a function of radius. Top Right: results for PL emissivity as a function of radius 
with a = — 1 (corresponding to q = —1). Bottom Left: results for PL emissivity as a function of 
radius with a = 0 (corresponding to g = 0). Bottom Right: results for PL emissivity as a function 
of radius with a = 1 (corresponding to g = 1). 


4.5. Correlation between L' and 7 ' 


Here we explore the implications of a possible correlation between F' and the random Lorentz 
factor of the electrons in the bulk frame 7 ' such that F' = Correlations of this kind may 

result in case the reconnection layer is continuously fed by incoming particles (i.e. k = 1) and 
a quasi steady state is achieved between incoming and outgoing particles. This behaviour has 
previously been found in PIC simulations (e.g., Cerutti et al. 2012[ 2013| ). There is a physical 
motivation for this since the electrons that are accelerated to higher energies are preferentially 
those that spend a longer time being accelerated in the reconnection layer, and thus their velocities 
tend to be more tightly collimated, down to a scatter of A6 ~ l/7e- This can be understood in the 
following way. The electric field in the reconnection layers is perpendicular to both the original 
magnetic field’s direction and to the incoming particles’ velocities. Electrons are only accelerated 
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as long as they have not defleeted signifieantly from the X-point and move along the eleetrie field 
E. If this aeeeleration is over a distanee l\\, then they aequire an energy of 7 'meC^ ~ eEl\i 
(3ineBol\i where Bq is the magnetie field strength in the bulk, before the reeonneetion layer. On 
the other hand, they exit the aeeeleration site after traversing a distanee l± in the perpendieular 
direetion, at whieh point their angle relative to the eleetrie field direetion is 


[7 rsj - - 

/|| 'y'^rrieC^ 




Rl{Bo) 7' 


(32) 


It follows that r' may be up to linear in 7 ' (see Eq. C.5 

0 < 7 < 1. 


in Q and that one might naively expeet 


In the presenee of sueh a eorrelation, eleetrons of different typieal energies may dominate 
the emission at different observed frequeney bands. This leads to a qualitatively different physieal 
seenario than that presented elsewhere in the text. The main differenee from the ease with no T'— 7 ' 
eorrelation has to do with determining the frame in whieh the eleetrie field vanishes. In ease 
a eorrelation between T' and 7 ' exists, eleetrons of different energies are experieneing different 
fields and, eontrary to the ease with no eorrelation, there is no single primed frame where the 
eleetrie field vanishes and where eleetrons ean gyrate freely around the magnetie field without their 
veloeities getting isotropized in the jet’s frame. Here the role of the vanishing eleetrie field frame 
is played by the jet’s frameQ As a result, in order for the eleetrons’ emission to remain beamed, 
they must radiate a signifieant fraetion of their energy before their veloeities get isotropized in the 
jet (rather than the plasmoids’) frame (i.e. within less than a single Larmour gyration). Therefore, 
sueh a eorrelation ean be expeeted to manifest itself only at very high frequeneies - at a band 
near the elassieal “maximal synehrotron energy”. We stress that the results obtained under these 
assumptions are therefore likely not to be relevant for understanding the observed properties of 
sub-MeV pulses. However, they may still have interesting observable implieations for the high 
energy emission. 


Notiee that the emitters with different 7 ' have different T' eorresponding to different eenter of 
momentum frames. Moreover, the average eleetrons’ Lorentz faetor in the emitters’ frame, 7 e > is 
related to 7 ' through T': 7 ; r' 7 " (see ^for a derivation). It follows that: 


7 


/ 

e 



r 7 


r' 


K^- 


ne 


1-77 


(33) 


In addition, if we assume that the emission is dominated by synehrotron radiation (as expeeted in 
a magnetieally dominated emission region), then = 27^7e^’ where B' is the magnetie field 
strength in the jet’s bulk frame (whieh is expeeted to be roughly eonstant for partieles of different 


'in the acceleration region (or X point), where > B^, there is no frame where the electric held vanishes. 
However, the radiation occurs mostly outside of this region, in the jet’s bulk frame. 
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energy, as they are all assumed to radiate in a similar loeation). It ean be related to the magnetie 
field in the emitter’s frame througlj^i?" ~ B'V. It is useful to relate the observed frequeney to F': 




r<n(7') 


TeB' 

2'KmeC \ K J 


(34) 


Sinee the peak speetral emissivity is determined in the jet’s bulk frame. The speetrum in this frame 
is given by the regular synehrotron formulas. As detailed above, a eorrelation between 7 ' and F' 
implies that the eleetrie field vanishes only in the jet’s (rather than the plasmoids’) frame and as a 
result we should eonsider eleetrons in the “fast eooling” regime (i.e. that eool signifieantly in less 
than a dynamieal time), and so the partieles’ speetrum is: 


die "" d\\ii 


Ax 


/-p 


(35) 


The resulting speetrum is therefore: (x u as for the ease of synehrotron without boosted 

emission in the jet’s bulk frame. 

Although the speetrum does not ehange for F' S> 1, there will be an angle dependent eut-off 
to the speetrum in this ease. Notiee that eleetrons with different Lorentz faetors 7 ' are emitted 
by emitters with different F'. Emission from an emitter moving at F' ean be seen up to an angle 
^obs ~ where cos 6 *Q^g = n' ■ (3', this implies a eut-off of the observed speetrum that is 

determined by the observation angle: 

TeB' 

= -'.UP = 1/9*J = ■ (36) 


As mentioned above, in order to maintain the anisotropy, eleetrons are doing less than one gyration 
around the magnetie field. Therefore, this effeet will take plaee at frequeneies above the “regular” 
maximal synehrotron frequeney r^syn.max- 


Another effeet introdueed by a eorrelation of the type F' cx H relates to the typieal width of 
the pulses. So long as the latter are determined by F' (see §4.4 for details), the typieal width will 
be different at different frequeneies (whieh are dominated by emission from eleetrons at different 


energies). As seen by Eq. 34 this will lead to AT cx u We return to this point in (5.3 


As is the usual ease for synehrotron, we ean safely assume that eaeh frequeney is dominated 
by eleetrons of a single energy]^ As mentioned above, the only frame with a purely magnetie 


^This is different from the usual relation that is derived under the assumption that there is a single primed frame 
where the electric field vanishes. 

^This is true so long as the spectrum obtained in this way is softer than the synchrotron spectrum for single 
electrons: oc . 
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field in case of a F' — 7' correlation, is the jet’s frame. We therefore start with the peak spectral 
emissivity of individual electrons in the jet’s frame (P^/ g aririeC^B'/ (3e) oc B') and use the 
regular Lorentz transformation (as in Eq. ([^, but with the direction of the transformation reversed) 
to obtain the (isotropic) luminosity in the emitters’ frame: 


L''„= dn' 






P' 


(37) 


z/' ,e,max 


din 7, 


syn/ 


/ e,syn/ 


'e,synp,/ , 7 5 

V /e,syn/ 


where 7e syn is the Lorentz factor (in the jet’s frame) of electrons for which the synchrotron fre¬ 
quency is the observed frequency w. 


( u27imecK 


[v{y)iy(y,4,)eB') 


1 

\ 2 - 7 } 


(38) 


Note that in Eq. (37) we have expressed the spectral luminosity in the jet’s frame {L'^,) in terms 
of the electrons’ energy distribution (dNe/d'j'J and the peak luminosity of individual electrons 
WhcrC /K'.e.ras.. ~ zfw (T^syn) ~ ^e(7e,syn) iS the nUmbcr of ClCCtrOnS 


/e,syn//“j/',e,max din 7 ' V 'e,syn/ ^'ev/e,syn/ 

whose synchrotron frequency is close to the observed frequency, so that each of these electrons 
radiates at that frequency near its peak spectral power P^/ e,max its peak isotropic equivalent 
spectral luminosity g ^^^. The latter are related through E'^, g ^^^ ~ T'^P^, ^ since P^, = 
dE'/dp'di' and L'^, = AT:{dE'/dv'dt'dVt') = ATT{dPl,/dVt') ~ (47 r/Af 2')P^, and in our case 
dvr/Afl' 4r'^ ~ r'^ as because of relativistic beaming, in the jet’s frame, the contribution 
of the luminosity is significant only within a cone of opening angle ~ l/F' in which V ~ F'. 

- NAP 

-‘■'ey /e,syn/ 


)^:^',e,max/r'(7e,syn) or ^at in the emitters’ 


/ e,syn 


The basic result of Eq. (37) is that E"„ 
frame the spectral power of a single electron whose synchrotron frequency is close to the observed 
one, satisfies ^ K"/^e{i^,syn) ~ ^^',e,max/r'(7;syn)> wMch is the usual familiar Lorentz 

transformation for the received emission in the jet’s frame, which corresponds to a single propa¬ 
gating electron, e.g. that is part of a blob and so corresponds to fc = 0. The correction for a steady 
state in the jet’s frame, fc = 1, is the extra power of 1 — A'ile^syn) sin 9' cos 0 (which is the ratio 
of the received and emitted times for a point source such as an individual electron) that appears in 
Eq. ( [39| ), which accounts for the fact that in this case there is no distinction between the emitted 
and received times, so that the radiation is not received over a shorter time than it is emitted, as is 
the case for a blob. 


In order to obtain the observed flux, we plug Eq. (37) into Eq. ([TS]) and integrate over 7': 


F.{T) = 


2r„7iP', 


V ,e,max 


(47r)2p2 


’2(m+l) Pymax 


dy 


'2/mm 


m+1 


m + y 
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_ 1 _ m 
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dcj) 
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/e,syn 
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(39) 




-{1-/3' (7e,syn) sin 9' cos 0) 
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The resulting speetrum is shown in Fig. The time integrated speetrum, as well as the speetrum 

2-p 

elose to the time of the peak flux approaeh iyF,y oc i/ 2 as expeeted for non-boosted emission. 
Results are plotted for p = 2.5 whieh is eonsistent with expectations from theory. Note that 
although for acceleration in reconnection sites a large range of p-values is found, depending on the 
exact setup and initial magnetization ( [Sironi & Spitkovsky |2014t[Guo et al.|2014[|2015||WernerH 
al.|2014[ Kagan et al.|2015p down to values of the order of p ~ 1.5, this of course cannot hold up to 
very large energies as they would require a divergence of the total energy, and is also inconsistent 
with both afterglow observations and observations of the high-energy spectral slope of the Band 
function. At either early or late times, and frequencies above z/max, the flux will become steeper. 


due to the angle-dependent cut-off described in Eq. 36 For these, the emission is dominated by 


large angles from the line of sight, with 6*T'(z/) > 1, and one can relate the observed spectrum 


F^{6') to the one along the line of sight, F,^{6' = 0) = F^{0) (e.g., Granot et al. 2002 20051, 


F,{e) [V\e') 


F,(0) V^^'(O) 


3-fc+f 




/2r,(fc-3-f) 


oc u 


-r;(3-fc+f) 


(40) 


F,{e) oc F,{0)u- 


' 2 ^ oc u 


Fig. [^clearly shows that at early and late times, the flux approaches this slope. The transition from 
the slope to this steeper slope, is at Umax- 


4.6. Results for the entire light-curve 


Under the same assumptions as in § 4.3[ and for the same dependencies of the luminosity on 
radius we calculate the resulting light-curve. In principle, the distributions of luminosities, Forentz 
factors and typical radii of emission, are model dependent, and are not uniquely determined by the 
fact that there is anisotropic emission in the jet’s bulk frame. These can reflect the inner workings 
of the central engine (see ^for more details). Since we do not want to invoke a specific central en¬ 
gine model, these distributions were chosen phenomenologically, motivated by typical light-curve 
observations. We assume that individual shells (leading to different pulses) are characterized by 
different typical duration Tq, different ejection times Tej, different T' and different peak luminosi¬ 
ties (i.e. different Fq). We assumed ten spikes in the light-curve, corresponding to ten emitting 
shells. In order to calculate the corresponding light-curve from ten shells, we randomly draw the 
parameters of each shell out of the following probability distributions. We draw the values of Tq 
for different shells from a log-normal distribution with a central value denoted as (Tq) and a typ¬ 
ical width of 0.3 dex. The ejection times Tej are randomly drawn from a uniform distribution in 
the range 0 < Tej < 20(To). We draw the values of T' for different shells from either a uniform 
distribution between 2 and 10 (or between 1 and 10). Finally, we draw the luminosities (T''„) from 
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Fig. 8.— Spectra of a single pulse at different times: T = [1.5, 2, 2.5]To, alongside the time- 
averaged spectrum, assuming a correlation of the type F' oc 7'’’. The electrons’ energy distribution 
is assumed to be ^ cc with p = 2.5. We also assume m = 0, fc = 1 and a log-normal 

emissivity with ax^R = 0.01. The peak of the light-curve is at T = 2Tq. The spectra is plotted 
for 7 = 0.8 (solid lines) and 7 = 0.5 (dashed lines). For comparison we also plot the expected 

2—p 

spectrum in this case, for non-boosted emission (and for the time-averaged spectrum): vFy cx z/~, 
and the asymptotic spectrum expected in this case for v > z/max- The frequency is in units of z/q, 
which is the typical frequency at which electrons with F' = 100 radiate (see 34), and below which 
the time-scale is dominated by the radial instead of the angular one. uFy is given in the units of 
z/qFo, where Fq is the peak of the time-resolved flux at z/q. 


a log-normal distribution with a central value denoted as {T"„) and a typical width of 0.5 dex. 
The results are plotted in Fig. For larger AR/Rq and for a = 0, the pulses tend to be wider 
(and the peak flux is weaker in order to conserve the total luminosity of the pulse). In this case, 
it is harder to resolve individual pulses and there is more overlapping of pulses. Generally, the 
resulting light-curves can be quite complex and have a wide range of typical behaviours, as seen in 
observations. However, it is basically a reflection of the choice of distributions above and therefore 
a true comparison with observations will depend on the specific model assumed. 


28 
















Fig. 9.— Light-curves of ten eonseeutive pulses with different L{R) (but the same for different 
pulses within a light-eurve). We assume a Band funetion speetrum and m = 0,k = 0. The 
flux is in units of the “typieal” peak flux Fo(r' = 1) as a funetion of the average normalized 
time T / (Tq). Tq, L"„ for different shells are log-normally distributed with widths of 0.3, 0.5 dex 
correspondingly. The ejeetion times are uniformly distributed in the range 0 < Tej < 20(To). 
Finally, F' is uniformly distributed in the range 2-10 (1-10) in the top (bottom) panels. Left: 
Power-law L{R) with different power-law slopes a and Rf = 2 Rq. Right: log-normal L(R) with 
different standard deviations ai^R. All parameters (except for F') are held constant for different 
light-eurves within any panel. The values of F' in the top (bottom) panel, from left to right, are: 
3.4, 9, 9.6, 2.6, 8, 7.6, 9.8, 5.5, 9.7, 6.4 (8,1.1,4.6, 7.4, 9.3,5.5,1.3,4.6, 6.7, 8.1). 


5. comparison with observations 
5.1. peak Luminosity-variahility correlation 


Several studies have explored the possibility that the variability of GRB light-eurves is related 
to their peak luminosities (Stem et al. 1999[ Fenimore & Ramirez-Ruiz 2000; Reiehart et al. 


20011. We use the definition for the variability presented in ([Kobayashi et al. ||2002|, whieh is a 


simplified version of that in (Reiehart et al. |200T I: 


EjFAT) - {F^{T))f 

Y^Fuiry 


(41) 
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where {F^{T)) is the flux smoothed with a boxear window with a time seale equal to the smallest 
fraetion of the burst time history that eontains a fraetion 0.45 of the total flux. [Reiehart et al. |([200T]) 


have found that the peak luminosities L are eorrelated with the variabi 

ity measures: L, 

Tiax 

y 3 . 3 lo 9 should however be noted that other studies JGuidorzi et al. 

2005 

Guidorzi 

2005 


have elaimed a eorrelation with a mueh shallower slope, of the order of l.SOlg 44 and 0.85 ± 0.02 


respeetively. The differenee between the results has been attributed to their respeetive treatment 
of sample varianee ( Guidorzi et al.||2006 ). These authors use a Monte Carlo method to show that 
when the sample varianee is of order the total varianee the latter methods beeome more reliable 
estimators of the slope in the Lmax — V eorrelation. 



log(V) log{V) 


Fig. 10.— Left: Maximal luminosity of a simulated light-eurve plotted against the variabilities 


assoeiated with the same light-eurves using a simplified version of the definition in (Reiehart et al. 


2001) assuming AR = 0.01i?o for ^ general emission meehanism (r) = 0) and for synehrotron 


emission with a eorrelation between T' and 7 ' as deseribed in § 4.5 The inner error bars on eaeh 


plot are the statistieal errors. Notiee, that these are small eompared with the intrinsie errors, as 
deseribed in §5.1[ In the latter, we explore different possibilities for the power-law index of the 
eorrelation, 7 . Right: Maximal luminosity as a funetion of variability for different assumptions 
on the width of the shells, T' (either varying independently between shells or eonstant within a 
burst), and the values of Fq for different pulses (varying log-normally or eorrelated with T'). 


A eorrelation between peak luminosity and variability may also be expeeted in our model, 
sinee the peak luminosity inereases with T' and the width of pulses deereases with T' so long as 
AR/R < l/T'. Notiee, however, that this trend is reversed for T' < 2 (see Fig. whieh may 
eause a negative eorrelation if the typieal values of F' are in the range: 1 < F' < 2. In Fig. 10 we 


plot the results for the Lmax - V eorrelation for simulated light-eurves in our model. In ealeulating 
the light-eurve we assume a typieal Band funetion speetrum emissivity (an = 0, /3 b = —2.3) as 
well as m = 0, fc = 0 and a Gaussian dependenee of luminosity on the logarithm of the radius. 


30 





































































We first let Tej, F' vary between different pulses (aeeording to the distributions assumed in 
4.6 see Fig. [^. Next, we eonsider synehrotron emission with a eorrelation between F' and 7 ' 


as deseribed in § 4.5 As mentioned above, sueh a eorrelation is less likely to hold for eleetrons 


radiating at frequeneies eorresponding to the sub-MeV peak. However, it still remains an open 
question whether a similar relation between variability and luminosity also exists at the GeV range. 
It is therefore important to provide theoretieal predietions for this situation whieh eould be tested 
by observations at a later stage. Again, we let Tej vary randomly for different pulses. We ealeulate 


the light-eurve at a eonstant observed frequeney. As ean be seen from Eq. 34 this means that we 
should eonsider 7 ' = constant between different pulses, implying K (xV. We therefore vary K 
using the same distribution of F' as for the ease with no eorrelation. In addition, we assume k = 1 
for this seenario, as deseribed in § 4.5 The light-eurve is then ealeulated using Eq. (39). In both 
oases, we find that for a\^R = 0.01 there is evidenoe for a positive eorrelation between Lmax and V 
with a power law index in the range 0.8—4.8 (see left panel of Fig. [T^, eonsistent with the observed 
values. To quantify the degree of eorrelation, we ealeulate the Pearson eorrelation ooeffioient for 
the generated light-eurve distributions deseribed above, and find it to be 0.7-0.9. The errors in 
Fig. 10 denote the smallest range within whieh lie 68.3% of the values of logio(Lmax) obtained 
from different realizations of the light-eurve. Note that the large errors in Fig. are systematie 


and not due to insuffieient statisties (as is shown explieitly in the left panel of this figure, where 
we eompare the systematie and total errors). They refleet the faet that the variability is set by a 
eombination of many different faetors, sueh as the exaet F' distribution for the different pulses and 
whether or not there are overlapping pulses in the light-eurve. For this reason a given degree of 
variability ean be obtained by light-eurves with different intrinsie properties, and different values 
of Lmax- Note also that the obtained distribution of Lmax does not follow a simple log-normal 
distribution and the errors in general may be asymmetrie. 

Variation in the dissipation rate and in the assoeiated eo-moving luminosity over AR <C 
R as refleeted in ainR = 0.01 may be present and have some eontribution to the overall light- 
eurve. However, it is probably more realistie to eonsider wider shells, or larger AR/R. When 


eonsidering a\^R = 0.1 in the right panel of Fig. 10 (and also by letting Fq vary log-normally. 


as in § |4.6[ and independently of F'), and leaving the other parameters as before, we find that the 
positive eorrelation is wiped out. This is beeause shells with the largest values of F' are no longer 
neeessarily the narrowest. In addition, whereas Lmax depends on the maximal F' in the burst, the 
variability depends on the width of all the pulses in the light-eurve. Sinee we have assumed no 
eorrelation between F' values in different pulses, it is natural that no strong eorrelation is seen 
between 

Nonetheless, reeonneetion simulations suggest that F' inereases with the magnetization pa¬ 
rameter a ( Guo et aL]|2014 , Werner et al.|2014 ), albeit very slowly, and it is not unreasonable that 
the typieal magnetization would vary between different GRBs. Thus motivated, we also explore a 


simple model in whieh F' varies between different GRBs as in § 4.6 (i.e. uniformly between 2 and 
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10) but is constant for a given GRB. We also explore the possibility that Fq ocT'. The eorrelations 
for these simulated light-eurves ean be seen in the right panel of Fig. 10 Even with these modifiea- 
tions, a positive eorrelation between variability and Lmax, for shells of width AR k, 0.1i?o is only 
obtained when Fq cx F' (resulting in a eorrelation with a power law index in the range 0.4 — 1.2 
and a Pearson eoeffieient of 0.3-0.6). This power law index is smaller than the observed value of 

3.3) ')gg reported by Reiehart et al. (2001), however it is eonsistent with the later measurements of 

1.3) ')o; 4 and 0.8 ± 0.02 found by (??) (see diseussion above). 


5.2. peak luminosity - peak energy correlation 


There have been various elaims in the literature that Up may be eorrelated with the peak lu¬ 
minosity, Lp, ( Yonetoku et al.||2004 ; Ghirlanda et al. ||2005 ; Yonetoku et al.||2010 |. The reported 
eorrelations vary in different studies and range between: Lp oc We aim to test the viability 

of sueh a eorrelation in the eontext of our model. 


First, reeall that: 


r'r. 


i VpLy^ 


iTFRVL') ^ V^L' 


(42) 


The faetor of F^ for the luminosity ean be understood as follows. The emitted power in the jet’s 
bulk and observer frames is equal, dE/dt = dE'/dt', but the reeeived power, dE/dT (where T is 
the observed time), is larger by a faetor of dt/dT = 1/(1 — f3p) ~ F^ in the observer frame. This 
eonsideration aeeounts for the angular distribution of the emitted radiation through the estimate 
dt/dT ~ F^, whieh effeetively assumes that the observed power is dominated by material whose 
radiation is beamed towards the observer. Alternatively, the emitted power per solid angle per 
partiele (Liso,e = dE/dtdVt) seales as ~ F"^, but the fraetion of the spherieal emitting shell size 
that is within the visible region is ~ 1/F^ so that together we onee more obtain a faetor of ~ F^. 
Combining the two relations in Eq. ([42]), we obtain: 




LL 






/2 


(43) 


This implies that if the ratio between the total luminosity in and square of the peak energy in the 
jet’s bulk frame, L'Jvp, is roughly eonstant, then a similar relation holds also in the observed 


frame, Lp oc al, whieh reproduees the observed eorrelation (see also Wijers & Galama 1999 


Ghirlanda et al. 2012). This ean naturally oeeur if the properties of the emission are set in the 
jet’s bulk frame. Even more importantly, the result is independent of F whieh is likely to vary 
signifieantly from burst to burst (and even within a single burst) and is unlikely to be direetly 
related to the properties of the flow in the eo-moving frame. 
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In case of anisotropic emission in the jet’s bulk frame, where the emission properties are 
set in the emitters’ frame, as eonsidered in this paper, the situation is more eomplex. It is most 
instruetive to examine the behaviour for F' S> 1. From Eq. (19 
l/F' (and motivated the definition of Fq cx 

this regime ~ FoF'^”^^ max(l, i?o/Ai?) where we identify 

T" 


that holds for AR/Rq > 


in Eq. |26|) it ean be seen that in 

Vp/TV and 




L ‘ 


L''lv'' On the other hand, Eq. (23) that is for a delta function emission in radius 


and reasonably holds for AR/Rq < l/F' shows that in this regime Ly^/L'l 
LpjL'^ ~ {v”lvp)L^jL'l„ ~ (Lj,^/L"„)/(FF'), this implies that 


FnF 


/2-2k 


Sinee 


r f T" 

_ ZA = P -p/-2fc J 

Z /2 Z /'2 Z /"2 ■'* ’ 


h 



AR/Rq > 1 , 

1/r < AR/Rq < 1 
AR/Rq < 1/r . 


(44) 


That is, L'//v'/‘^ = Lp/u^ = 1 is satisfied (i.e. this sealing extends to the emitters’ 

frame) only for k = 0 (blob-like emission) and AR > Rq, but not for k = 1 (eontinuous emission) 
or AR < Rq. In the latter eases a residual dependenee remains on F' (for AR < Rq), and 
sometimes also on AR/ Rq (for l/F' < AR/Rq < 1). Nevertheless, it is still reasonable that some 
of the properties of the emission in the emitters’ frame would be determined by proeesses that 
oeeur in the jet’s bulk frame, in whieh the emitters are energized and launehed at F', so that they 
eould in prineiple eorrelate with F' and still produee the observed eorrelations. 

We stress that whether or not the ratio L'p/^p is indeed eonstant, depends on the specifies 
of the emission meehanism, whieh are not uniquely determined by the model at hand, namely, 
anisotropie relativistic velocity distribution of particles in the jet’s bulk frame. It is nonetheless 
intriguing that if such a correlation indeed holds, then it will not be wiped out by variability in F 
(and possibly also in F'). 


5.3. spectral evolution in a single pulse 


Using the results of 


4.2 we can estimate the typieal width of a spike in the light-eurve at 
different frequencies. As ean be seen in Eq. (22), so long as there is no 7 ' — F' eorrelation, the 


width of the pulse is mostly independent of frequency. However, switehing between frequeneies at 
different parts of the spectrum (say from the peak to the lower energy speetral slope) ean somewhat 
ehange the pulse width. 


In Fig. 11 we plot the resulting light-eurve at different frequeneies, for a Band speetrum 
emission. We assume here a typieal Band funetion (a^ = —l,/3s = —2.3), sinee the intrinsie 
speetra is model dependent and not uniquely defined by the anisotropie emission model. We note 
however, that the results presented in this section are not strongly dependent on this assumption. 
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Fig. 11.— light-curves of a single pulse at different frequencies. All cases are plotted with a 
typical Band function {as = —^,Pb = —2.3), m = 0 and a log-normal dependence of luminosity 
on radius with AR/Rq = 0.01. The maximal flux is normalized to 1, and the time is in units of 
the typical time Tq and starting at Tq (when the first photons reach the observer). Left: No 7 ' — F' 
correlation. In this case we assume F' = 3, fc = 0 and a typical Band function. The shape of 
the peak is roughly constant at different frequencies, but tends to becomes slightly broader below 
the peak, right: A 7' — F' correlation with 7 = 0.8, k = 1 and with the “typical” values of F' 
corresponding to the observed frequency (see Eq. [^, being: 3.4, 8.5, 21.5, 54 from the lowest to 
highest frequency. In this case, the spectrum is (see H.S). Here the pulse clearly becomes 
narrower with lower frequencies, approximately following the relation AT cx implied by 
Eqg 


since they are mainly governed by the peak frequency and peak flux which change more as a 


function of time (as will be shown in Fig. 13). The light-curves of a single pulse at four different 
frequencies compared with the peak: u = O.Olz/p, O.lz/p, Up, lOup. All cases are plotted with F' = 
3, m = 0 and a Gaussian dependence of luminosity on the logarithm of the radius ai^R = 0.01. 
For the case with no F' — 7' correlation, the shape of the peak is roughly constant at different 
frequencies, but tends to becomes smoother and somewhat broader at lower frequeneies. However, 
if a correlation between F' and 7' exists, then higher 7' can correspond to a larger u^yn and a larger 


F' and thus a narrower spike, of width (see Eq. 

3^ 

i) AT cx l/F' cx 

power law index of ~ 0.4 — 0.5 ( 

Eenimore et al. 

1995 

Norris et al. 


would suggest 7 ~ 0.8 — 1. However, we stress that due to isotropization of the electrons’ orbit 
due to Larmour gyration, this effect is expected to be mostly dominant near z^syn.max ~ few x GeV, 
whereas in observations, the pulse widening can be seen at much lower frequencies, down to X- 
rays. 


Another result of this is a change in the spectral shape over time. In Fig. 12 we plot the spectra 
of a single pulse at nine different times: T = [1.01,1.1,1.5, 2, 2.5, 3, 5, 7, 9]To for m = 0, fc = 0, 
a constant dependence of luminosity on radius from Rq to Rj = 2 Rq, and for both F' = 1 and 


34 




































spectrum at different times, F' = 1 




Fig. 12.— Spectra of a single pulse at nine different times: T = [1.01,1.1,1.5, 2,2.5, 3, 5, 7, 9]To 
for r' = 1 (left) and F' = 3 (right). The intrinsic spectra is assumed to be a typical Band function 
(ttB = —1, (3b = —2.3) with m = 0, k = 0 and a constant emissivity from Rq to Rj = 2 Rq. The 
peak of the light-curve is at T = 2To. The spectra is plotted for both T' = 1 (left) and T' = 3 
{right). uFy is given in the “typical” units of UpFo and the frequency is in units of the intrinsic peak 
frequency Up. The spectrum evolves from soft and weak to hard and strong and back again. 


T' = 3. The spectrum evolves from soft and weak to hard and strong and back again. Observations 
of GRB pulses have shown two typical behaviors: intensity tracking and hard to soft evolution 
( |Ford et aL~||1995[ [Preece et al. ||2000] [Kaneko et al. ||2006 ; |Lu et al. ]|2012| ). Indeed, Figure 
(as well as Figs. [T3] and [14]) demonstrates that our model predicts intensity tracking evolution for 
sufficiently large F' (approximately F' > 2), while a hard to soft evolution appears for low enough 
r' (approximately F' < 2). Since pulses in our model are asymmetric, with the rise time much 
faster than the decay time, this will make the hard to soft trend much easier to detect than the initial 
soft to hard (or intensity tracking) evolution, and may account for the observed population of GRB 
pulses exhibiting “hard to soft” evolution (for instance this was shown for GRB 130427A in the 
best analysed pulse to date ( Preece et al. [2014 1). Note that a hard to soft evolution is also obtained 
for F' < 2. The evolution of the spectral parameters can be seen more directly in Fig. 13 where 
we plot Up and the Band photon indices aB,(3B as a function of time. The latter do not change 
significantly in time and remain constant on their intrinsic values. We also note that the change in 
spectrum becomes significantly more pronounced for as F' increases. For the latter, the softening 
and weakening of the spectrum in time follows the behaviour predicted for regular high-latitude 
emission: F^, = where a is the relevant spectral slope of F^, depending on the observed 

portion of the Band spectrum (either a = — (a^ -f 1) or a = —(/3b + 1)). This behaviour can 
be seen more clearly in Fig. 14 for F' = 1 and to a lesser extent also for F' = 3, 10 at late times 
(which would be harder to observe, as the flux is by then much lower than its peak value). Initially, 
the decay for F' = 3 is much steeper and might be able to account for the observed steep decay 
phase at the beginning of GRB afterglows (e.g., |Tagliaferri |2()05| ) . 


35 








































Fig. 13.— Spectral parameters of a single pulse at different times (Left: T' = 1, Center: F' = 3, 
Right: r' = 10). The figures are plotted with a constant emissivity from Rq to Rf = 2 Rq. The 
peak frequency is plotted in red and the Band photon indices as, [3 b are in black (solid and dashed 
lines accordingly). The intrinsic spectra is assumed to be a typical Band function (as = —1, [3b = 
—2.3) and m = 0, /c = 0. At first the peak energy increases for a short while, and later it decreases 
significantly. This behaviour is more pronounced for larger values of F'. The spectral indexes are 
basically equal to their intrinsic values and are almost constant in time. 


6. Discussion 


In this work we have considered a model for a relativistic source that emits anisotropic ally in 
its own bulk rest frame. Our main motivation was to calculate the most prominent expected obser¬ 
vational signatures for prompt GRB emission in which the dissipation that powers this emission is 
driven by magnetic reconnection. This motivated our choice of the geometry of the magnetic re¬ 
connection sites (thin spherical shells - reconnection layers) and the bulk Lorentz factor (F ^ F'). 
This model may, however, also be applied to other anisotropic emission models in GRBs, as sug¬ 
gested for instance in the relativistic turbulence models and the “jet in jet” models. This model 
may also be relevant to other astrophysical systems, such as active galactic nuclei, where very 
short time-scale variability (minutes) has been discovered (e.g., Aharonian et al.||2007 ; band et al. 


20071 possibly suggesting emission from mini-jets within the main relativistic jet (Giannios et al. 


20091. Another possible context for relativistic motions in the bulk frame are the GeV flares in the 


Crab nebula. These flares were suggested to be Doppler boosted emission from the termination 
shock ( {Komissarov & Lyutikov |2011| ) or on the hot spot at the base of the jet feature in the nebulae 
( {Lyubarsky 2012| ). In both cases, however, the bulk frame is not moving relativistically towards 
us, but is instead effectively at rest. 


We have considered a thin spherical shell propagating at a Lorentz factor F 1, which emits 
as it expands outwards and either (i) turns on and off abruptly at radii Rq and Rj = Rq + AR, 
respectively, with the co-moving luminosity scaling as a power-law (of index a) with radius within 
this range, or (ii) turns on and off more gradually where the co-moving luminosity has a log- 
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Fig. 14.— peak of z/Fj, (normalized sueh that its maximum value is 1) as a funetion of Up (normal¬ 
ized sueh that the initial peak of the speetrum in the bulk frame is 1) for a single pulse at different 
times (F' = 1 - solid line, F' = 3 - dashed line and F' = 10 - dashed-dotted line). The figures are 
plotted with a typieal Band funetion (as = —1, = —2.3), m = 0, k = 0 and a eonstant emis- 

sivity from Rq to Rj = 2 Rq. Arrows denote the direetion of progressing time along the eurves. 
r' = 1 shows a “hard to soft” behaviour, whereas F' = 3,10 show an “intensity traeking behaviour. 


normal dependenee on radius (eentred on Rq with a standard deviation ainij;). The emitting material 
is assumed to move at a Lorentz faetor F' relative to the shell’s loeal rest frame, in a direetion 
transverse to the shell’s veloeity (i.e. normal to the radial direetion). This direetion of motion is 
assoeiated with the anti-parallel reeonneeting magnetie field lines, and equal amounts of material 
are assumed to move along the eorresponding two opposite direetions (e.g., on either side of an 
x-point within the thin reeonneetion layer). 


Strong magnetie fields may also play an important role in the dynamies of relativistie jets. 
In AGN jets, the Thompson optieal depth at their base is insuffieient for effieient thermal aeeel- 
eration by radiation pressure, suggesting that Blazars’ jets must be magnetieally aeeelerated and 
therefore initially Poynting-flux dominated. In GRBs thermal aeeeleration is a viable alternative, 
but hydromagnetie jet launehing is expeeted to dominate, and the required high magnetization near 
the souree helps avoid exeessive baryon loading, whieh might prevent the jet from aehieving the 
required high-Lorentz faetors inferred from the GRB prompt emission (for a reeent review on the 
role of magnetie fields in GRBs see Granot et al. 20151. Uneonfined steady-state, ideal MHD 
axi-symmetrie outflows that start highly magnetized (with (Tq 3> 1) near the eentral souree remain 
highly magnetized (cr ~ far from the souree, while eonverting only a small fraetion of their 
initial eleetromagnetie energy to bulk kinetie energy (F ~ ). This may be alleviated if the 

outflow is gradually eollimated into a narrow jet of half-opening angle 9j, whieh eould asymptot- 
ieally lead to a ~ and F ~ or due to an abrupt drop in the external pressure 
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(e.g. when a long-duration GRB jet breaks out of its progenitor star), but in all eases the asymp- 
totie magnetization is still fairly high, a > 1. However, the strong Inverse Compton eomponent 
observed in Blazar speetra suggest a low a in the emission region (see e.g.|Ghisellini & Taveeehio 


2009 1), while in GRBs different lines of evidenee suggest a similar eonelusion (e.g.,|Beniamini &: 


Piran 20141. This result is eommonly refereed to as the a problem (Goldreieh & Julian 19691. 
This problem may be solved by relaxing one or more of the underlying assumptions of steady- 
state, ideal MHD and axi-symmetry. Strong time-dependenee ean lead to a mueh more effieient 
eonversion of magnetie to kinetie energy, and to a < 1 far from the souree that allows effieient 
dissipation in internal shoeks ( Granot et al.||2011[ ). Alternatively, magnetie reeonneetion ean play 
an important role in the jet’s aeeeleration and in the reduetion of a, and may be greatly enhaneed 
by various instabilities, some of whieh or non-axisymmetrie (sueh as the kink instability). The dis¬ 
sipation of magnetie energy generates heat, some of whieh ean be ehannelled to radiation, but most 
of whieh ean eontribute to the (now thermal) aeeeleration of the jet. The distanee from the eentral 
engine where effieient reeonneetion takes plaee, will determine whether the jet is still aeeelerating 
([Drenkhahn fc SpruIt~||2002[|Lyubarsk^|2010[), eoasting, or deeelerating (|Lyutikov & Blandford 


2003 [ IGiannios ]|2006| ) at that stage. In GRBs sueh late reeonneetion during the jet’s deeeleration 
phase might be expeeted for the following reason. The jet loses lateral eausal eontaet during its 
aeeeleration, eventually reaehing T9j ^ 1, so that onee it starts being deeelerated by the external 
medium inereasingly larger angular seales (~ l/T) in the jet eome into eausal eontaet and allow 
reeonneetion of magnetie struetures on sueh seales, whieh were not possible at earlier times. 

Effieient magnetie reeonneetion may naturally be aehieved in a striped wind magnetie field 
geometry. In this ease, the light-eurve variability time may refieet the time for the flipping of the 
magnetie field at the souree (i.e. the base of the outflow) over the (dimensionless) eo-moving flow 
veloeity into the reeonneetion layer, For a magnetar eentral engine, a periodie flipping of 
the field’s direetion is expeeted with a sign flipping time equal to half of the magnetar’s rotational 
period P, as the loeal field switehes sign twiee in eaeh rotation (e.g., Usov |1992). This may lead to 


relatively “ordered” light-eurves with short variability time-seales, AT 5(P/1 ms)(/9in/0.1) ^ ms. 
These beeome even shorter for a relativistie inflow veloeity ~ 1 into the reeonneetion layer that 
may arise for large cr-values (e.g., |Guo et all]|2015| ), whieh also leads to aeeeleration of the emit¬ 
ters to larger values of T'. Sueh short variability times are typieally signifieantly shorter than the 
dominant variability time in most GRB light-eurves. 

An alternative option that ean naturally lead to larger variability times as observed in most 
GRB light-eurves is that the eentral engine is an aeereting blaek hole, where the magnetie field that 
is adveeted inwards with the aeereted material randomly switehes sign either through some disk 
instability or by aeereting blobs of plasma with a randomly oriented frozen-in magnetie field (e.g., 
Blandford fc Payne~||1982[ |Rees |[T9Mt |Uzdensk)^ fc MaeFadyen |[2006 [Barkov fc Komissarov 


20081. In this ease the typieal time-seale for a flip in the magnetie field orientation is more uneertain 


and in addition, due to the stoehastie nature of the proeess, the resulting seale for the flipping of 
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the magnetic field will be much more variable. Therefore, it can potentially account both for the 
much larger observed typical variability times compared to the naive expectation for a millisecond 
magnetar, as well as for the large range of variability time-scales that are observed in most GRB 
light-curves and their stochastic nature in general. The resulting efficiency of reconnection in 
such an outflow with a random flipping of the magnetic field is uncertain, and may be somewhat 
lower than for a periodic flipping that was mentioned above. This might make it hard to satisfy 
the high gamma-ray emission efficiencies that are commonly inferred from GRB observations 
( [Panaitescu and Kumar]|2002| [Granot et al.||2006] [Fan et al^|2006| ). However, as was recently 
shown by Beniamini et al. ( 2015| ) these efficiencies should actually be only moderate (~ 0.1), 
significantly less than what was suggested by some earlier estimates (~ 0.9). 

Throughout most of this work we have focused on the properties of the emission that are 
direct results of anisotropic emission in the jet’s bulk frame, and that are largely independent of 
the specific radiation mechanism/s at work. However, as mentioned above, our main motivation 
comes from magnetic reconnection models for the energy dissipation. In reconnection models 
for the prompt emission of GRBs, because of the large expected magnetic fields the most natural 
radiation mechanism is synchrotron. We consider the emitting electrons to be in the “slow cooling” 
regime in order for the spectrum to be consistent with observations in the optical and X-ray bands 
( Beniamini fc Pirarr||2014 ). This can be expected if the frame in which the electromagnetic field 
becomes purely magnetic, is that of the emitters (“the plasmoids”). However, if the latter is in fact 
the jet’s frame, then electrons producing the observed sub-MeV peak will isotropize much faster 
than they can radiate their energy and the resulting emission would no longer be anisotropic, except 
at frequencies very close to z/gyn^max (see ^4.5| for details). In “one zone” magnetically dominated 
emission regions, the emitting electrons would not be in the “slow cooling” regime, unless T > 600 
( Beniamini k. Pirarr||2014 |; and even then, slow cooling would only occur for the largest possible 
radii, and mildly relativistic electrons. However, in the anisotropic model, suggested in this paper, 
there are two important differences compared to the “one zone” case. First, as mentioned in ^ 
the radius implied by the variability time-scale is larger by a factor F' as compared to the non- 
boosted case (leading to a weaker magnetic field for the same overall luminosity). Second, the 
local magnetic energy density within the reconnection sheet, where the particles are emitting, can 
be up to an order of magnitude lower than the average field. Combining these two facts implies 
that in these models, even for a smaller bulk Lorentz factor: F 300, electrons with 7 ' < 10 will 
be in the ‘slow cooling” regime, as required by both theory and observations. In addition the upper 
limit on 7' will increase significantly (7' oc F®) for larger values of F. 


We briefly summarize our main results and their comparison to the observations: 


1. Light-curve Variability : GRB light-curves show highly variable time profiles with the 


times between pulses roughly equal to the pulses durations ( |Nakar & Piran |[2QQ2| ). We 
show that isotropic reconnection emission models are in conflict with this observation since 
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they generally imply pulse widths signifieantly broader than the time between pulse peaks. 
However, we find that the anisotropie emission model presented in this work ean produee a 
highly variable emission at large distanees from the eentral souree, as required by observa¬ 
tions. 


2. Pulse Asymmetry : observed GRB pulses tend to be asymmetrie with a typieal rise to deeay 


time ratios frise/^decay ~ 0.3 - 0.5 ( 

Nemiroff et al. ||1994 

Fishman & Meegan||1995 

Norris 

et al. 1996[ Quilligan et al. 2002 Hakkila & Preeee 

20111). Asymmetrie pulses with 

frise < fdecay naturally oeeur in our model (see 1 

4.4 

and Table with f rise/f decay depending 


on the radial dependenee of the bulk Lorentz faetor (F oc and eo-moving luminosity 

(power-law oc or log-normal f(R/Ro)), on F', and on the radial extent of the emission 
produeing the pulse, AR/Rq. In partieular, they naturally oeeur for AR > Rq (e.g. when 
f{R/Ro) is asymmetrie, see for instanee the log-normal or power-law funetions with a < 0, 
eonsidered in this paper). Sinee the shape of a pulse in this ease is determined by f{R/Ro), 
observations of the pulse strueture eould be used to deduee the evolution of the emissivity 
powered by magnetie reeonneetion at a given site in the jet, in a way that may be largely 
independent of the miero-physieal properties of the reeonneetion and instead depend mainly 
on its gross properties. We note however, that sinee we effeetively see radiation from only a 
small angular portion of the jet, it would be hard to use this in order to learn about the global 
strueture of the jet. The pulse-shape gradually beeomes more symmetrie as F' inereases, so 
long as 1/F' > AR/Rq. At the same limit, the pulse both beeomes narrower and peaks 
at higher fluxes (for a given total radiated energy) for larger F'. Finally, we note that for 
F' > 2 and AR/Rq < 1, the degree of asymmetry depends mainly on AR/Rq. Here, onee 
more, observations of the pulse asymmetry eould be used to deduee the radial profile of the 
emission in the jet as mentioned above. 


3. Luminosity-Variability Correlation : previous studies have found a eorrelation between 
the variability of GRB light-eurves and their peak luminosities ( Stern etal. |1999 ; Fenimore 
I&: Ramirez-Ruiz ||2000[ Reiehart et al. ||^01t [Guidorzi et al.||2005 Guidorzi 2005 1. For 
AR/R < 1/F' and F' > 2 our model prediets that as F' inereases, the peak luminosities 
of pulses inerease while their widths deerease (leading to larger variabilities). For shells of 
width AR/Rq < 0.01 or even (assuming a natural assumption that the luminosities are eor- 
related with F') AR/Rq < 0.1, a eorrelation between maximal luminosity and the lighteurve 
variability is indeed obtained in our model. Although in the latter ease, the power law index 
of the eorrelation is signifieantly smaller, it is still eonsistent with the values reported by 
( [Guidorzi et al.|2005||Guidorzi |2005] ). 


4. Lp-z/p Correlation : many studies have elaimed the existenee of a eorrelation between the 
peak luminosities Lp and peak frequeneies Up of GRBs ( Yonetoku et al.|| 2004, Ghirlanda et 

Reeently, it was further suggested 


al. 


2005 


Yonetoku et al. 


2010) sueh that L 


, 1.5-2 

p i^p 


40 








































































( Guiriec et al.||2015 ) that these results eould hold also within a single burst, eomparing dif¬ 
ferent pulses. In both the non-boosted (F' = 1) and anisotropie models, both the peak 
luminosity and peak frequeney are Doppler boosted from the emitters’ frame. This results in 
a relation between the luminosity and the peak frequeney in the observer frame that is sim- 
ilar to that in the jet’s bulk frame: independent of the bulk Lorentz faetor F. This 

relation implies that if sueh a eorrelation between peak frequeney and the luminosity exists 
in the jet’s bulk frame, then it would naturally lead to a similar eorrelation in the observer 
frame, and would not be wiped out by varying F, whieh is likely to ehange signifieantly 
from burst to burst (and even within a single burst) and is unlikely to be direetly related to 
the properties of the flow in the jet’s bulk frame. In ease of a blob-like anisotropie emission 
with AR/Rq > 1/F' this result ean be extended further as it holds also in ease the emission 
properties are set in the emitters’ frame. 

5. Pulse Widths and Spectral Lags : observations of GRB light-eurves at different frequen- 
eies show that the pulse widths tend to deerease with frequeney approximately as IF(z/) oc 


u 


-0.4 


( Fenimore et al. |1995t Norris et ah 1996i 2005t Bhat|2012 ). A related observation is 


that at larger frequeneies, pulses tend to peak earlier - this is usually referred to as speetral 


lags (Norris et al. 1996; Band 19971. We note that in ease the pulse shape at different 


frequeneies is the same, and the emission starts at the same time, this would follow direetly 
from the anti-eorrelation between pulse width and frequency mentioned above. In our model, 
in case there is a correlation between the bulk Lorentz factors of the electrons,7', and F', then 
it is possible to obtain a negative correlation between the pulse width and the observed fre¬ 
quency, similar to the observed one. This is because higher 7' typically correspond to both a 
larger Usyn and a larger F' that leads a narrower spike. The observed relations are reproduced 
for F' oc 7'° ®“^. However, we stress that due to isotropization of the electrons’ orbit due to 
Larmour gyration, this effect is expected to be mostly dominant near i^syn.max ~ few x GeV, 
whereas in observations, the pulse widening can be seen at much lower frequencies, down 
to X-rays. 

6. Single Pulse L — i/j, Correlation : observations of GRB pulses have shown two typical be¬ 
haviours for their evolution in time: intensity tracking, where the peak frequency follows the 
overall flux or luminosity, and hard to soft evolution, where the peak frequency is constantly 
decreasing in time ( [Ford et aTr~||1995[ [Preece et al. ||2000| [Kaneko et al. ||2006t |Lu et al. 


20121. Interestingly, both behaviours are possible in our model, depending on F'. For small 


anisotropy (F' < 2), the the peak frequency always monotonically decreases, while the flux 
first increases and later decreases. Therefore, one may expect a “hard to soft” evolution in 
this case. However, for larger values of F' (> 2), the peak frequency also increases at first 
as the flux rises towards the peak of the pulse, and only later decreases during the tail of the 
pulse, so in this case one expects an “intensity tracking” behaviour. 
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7. Rapid Decay Phase : Finally, observations of GRBs in the early afterglow phase exhibit a 
so ealled “rapid deeay” phase ( [Tagliaferri [2005 [ ). This early phase of rapid flux deeay at the 
end of the prompt emission was unexpeeted sinee the flux often falls faster than would be 
expeeted from high-latitude emission (whieh is the shortest time-seale on whieh the flux ean 
deeay, even if the emission shuts off instantly). A possible solution in many cases is that the 
initial fast decay is dominated by the last pulse so that the zero-point for the power-law decay 
should be taken as the onset of that last pulse rather than of the whole GRB, i.e. the onset of 
the first pulse. In the anisotropic model the initial flux decay of a single pulse is significantly 
more rapid than in the non-boosted case (and increases with F'), and is thus able to account 
even for the most extreme cases of such rapid decay phases. A similar solution, producing 
rapid decays due to anisotropic emission was considered by |Beloborodov et aL~| ( |2011| ) who 
suggested an anisotropic emission model for GRB afterglows. 


We conclude that the simple anisotropic emission model presented in this paper is able to 
reproduce many of the observational features seen in the prompt emission of GRBs, and therefore 
seems very promising. In a future work we plan to extend this model to more realistic magnetic 
reconnection configurations, such as those obtained in relativistic reconnection simulations. This 
would help determine some of the free parameters that still exist in this model and could provide a 
self-consistent model for energy dissipation and emission in the prompt phase of GRBs. 


We thank Tsvi Piran, Daniel Kagan, Pawan Kumar and Wenbin Lu for helpful discussions. 
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supported in part by the ISF grant 719/14 (JG). 
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A. approximate expression for the flux 


Here we derive the result in Eq. ( [T9l ), by solving the double-integral in Eq. ( [T8| ) in the limit 
E' S> 1. Eirst, the inner integral, may be written as 

p27V 

h= (70(1 —/S'sin 0'003 0)^“““^ . (A.l) 

Jo 

Sinee the power k — a — 3 is typieally negative, the integrand beeomes maximal when (1 — 
/S'sin 0' cos 0) is minimal, i.e. for 0 = 0,6^' = 7r/2, and it rapidly drops for different values of 
(j), 9'. Therefore, in order to evaluate the inner integral we ean use the approximation (/> -C 1 
and coscj) ~ 1 — 0^/2, and take the limits of integration to be form — cx) to cxd (as the dominant 
eontribution is anyway from a small region where |0| < l/T' 1), 


h 



1 -/S'sin^' 



k—a—3 


A/^r(2.5-f a — fc) (1 —/7'sin 6*')^ “ 

T{3 + a — k) a//S' sin 9' 

(A.2) 


This preprint was prepared with the AAS lAT^ macros v5.2. 
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where here T is the gamma funetion. As for the outer integration, over y, sinee most of the eontri- 
bution is from 6' ~ 7r/2, this eorresponds to ~ j/o = (m + so that 


m+1 


2+a 


y 


a-l-f (l+«) 


1—2—a 


m(l+Q;)+ 2 (l —a.) 

[m + 2) 2 (m + l) 


(A.3) 


Also, the integral does not vanish only if the eorresponding radius, R 
2/oi?o(T/To)V(™+i) is in the emission region, i.e. between Rq and Rf, or 


1 < 


T 


_{m + 2)To 


1 

m + 1 
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Ro 


vRl ~ VoRl = 


(A.4) 


Moreover, -y/ [3’ sin 9' ~ 1, so that the only term in the integrand that signifieantly varies in the 
region where the dominant eontribution eomes from is (1 — (3' . Now we ean define 

y = y — yQ and approximate 


. + 1-1) 
sinb* = 2— -;- 1 


m + y 


—m—1 


8^2(m+2) • 


(A.5) 


Taking also the limit where T' 3> 1 so that ~ 1 — l/2r'^, we have 
{I - (3'^ 




2P2 


(A.6) 


Changing variables io y = yV/2y'^^‘^, we have dy = dy = {2y^~^‘^/r')dy. Sinee the main 
eontribution to the integral is from y of order unity, then for T' 3> 1 we ean take the limits of 
integration over y form —oo to oo, and when the eondition in Eq. (|A.4|) is satisfied one obtains 




dy{l — /S'sin 0 
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' J/min 
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(A.7) 


Putting all the terms together we obtain 
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B. Angular time-scale in the anisotropic emission model 

We aim to find here the angular time-seale for the model with anisotropie emission presented 
in this paper and show that it is Te = (1 + z)R/cV‘^V'. We define ^ = (r6*)^. In the relativistie 
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limit, r 3> 1, the main contribution to the radiation comes from 6 <^1 for which cos 6* ~ 1 — 0^/2, 
and /? = n/c ~ 1 — l/2r^ so that 


R= f3{t')cdt' ~ ct— 


cR 


dR' 


'o 212 (i?') 


= ct- 
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2(m + l)r2(i?) 
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t ^ — 
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1 + 


2{m + l)r2 

(B.l) 

Using the definition of the EATS along with the relations R/Rq = and R/T"^ = 

{R/Ro)'^^^Ro/Tl = y'^^^{T/To)Ro/Tl (see Eqs. [|^ and [ 101) one obtains: 


T = {l + z) 


Rcos9 


(1 + z)R 


[1 + (m + l)e'] = Ty^+^ [1 + (m + 1)^'] . (B.2) 


2(m + l)cr2 

This can be inverted to obtain an expression for ^ (see for instance ([Genet fc Granot |2009|)): 
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(B.3) 


Plugging this into the expression for sin 9' (Eq. [T^ we obtain: 

4(m +l)(?/-”^-i - 1) 4^ 


• 2 nl 

sm 9 = 


{m + y 


-m-l\2 


(1 + 0^ ■ 


(B.4) 


Eet us also define the angle between n' and jd' as x'- By definition, and using the fact that 
cos^Y = cos^Y', we have: 


cos^ Y = {n' ■ Y'Y = sin^ 9' cos^ 0' = cos^ 0 . 


(1 + 0^ 


Defining the parameter 

we thus obtain a quadratic equation for 

4/e = (1+0' 

for which the solutions are: 


/ = 


COS^ 0 
COs2 x' 


e' + O2-4/) + l = 0, 


(B.5) 


(B.6) 


(B.7) 


e ± = 2 / -1 ± 0 ( 2 /- 1 ) 2-1 . 


(B.8) 


The dominant contribution to the emission arises from within the beaming cone, i.e. from X' < 
l/V. Therefore, in order to find the boundaries of this region, we simply substitute x' = l/B' in 
the expression for the parameter / in Eq. (B. 61 . Hence, each value of / corresponds to a particular 
value of cos2 0, which in turn correspond to 4 different values of 0 = ±0o and tt ± 0o- For each 
of these 0 values, there are two values of ^ according to Eq. (|B.8), which each correspond 
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to a value of 0. As shown in Fig.[^ for a given value of F' this implies two approximately eireular 
regions of angular radius 1/(FT) around the two points of interseetion of the x'-axis (0 = 0, vr) 
with the eirele of angular radius 6 = 1/T around the line of sight (i.e. ^ = 1), i.e. (j)) = (1, 0) 
and (1, vr) or {6, cj)) = (l/F, 0) and (l/F, vr). For a given value of x' = 1/F', there are solutions 
for ^ only for / in the range !</<!/ cos^x^- The outer edges of this regions in the 0 
direetion eorrespond to / = 1 and (j)o = x' = l/F'. The eorresponding angular separation or 
distanee from the eenter of the relevant region is ~ 0o/r = l/(r'r). Along the x'-axis 0o = 0, 
/ = 1/cos^x^ = l/cos^(l/r') 1 + 1/r'^ or 2/ — 1 1 + 2/r'^, where the last approximation 

holds for r' 3> 1, for whieh in this ease one obtains 






for whieh the eorresponding values of 9 are: 

r0± = 
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i±- 
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f ^ FP 


9± ^ — ± 


(B.9) 


(B.IO) 


whieh again eorrespond to an angular separation or distanee from the eenter of the relevant region 
(that eorresponds to 0 = l/F with the same value of 0 = 0 or tt) of ~ 0o/F = l/(rT). 


We see that the typieal range of 9 from whieh a dominant eontribution to the emission arrives 
is indeed l/(rr'). We ean now estimate the angular time-seale, whieh is the differenee in arrival 
times of photons from 6 *i = l/F and 6*2 = l/F + l/(rr'): 


Te Rf ^ 

-= —(cos o'. 

1 + z c ^ 


D n2 T3 

oos9^)^^{^-^) = ^^{9^ + 9.){9. 


9.] 


2R 

cT'^r 


(B.ll) 


This is taking the diameter of this approximately eireular region. A better estimate is obtained by 
using its angular radius, whieh results a a time smaller by a faetor of 2, 


Te 


(1 + z)R 

cT^r 


(B.12) 


that is shorter by a faetor r '/2 than the angular time for the ease of non-boosted emission in the 
jet’s bulk frame. 


C. The typical electrons’ energy in the emitters’ frame 

Consider a distribution of eleetrons, sueh that for a given 7', eleetrons veloeities are distributed 
uniformly within a eone of half-opening angle 9o{'ye) around the x'-axis in the jet’s bulk frame. The 
four veloeity of an eleetron with a random Lorentz faetor 7' travelling at some angle 9 to the x'-axis 
in the jet’s bulk frame, is given by: 

= (7e, 7e^e/i, 76/56^1^1^, 0) (C.l) 
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where /3e = \/l — 7 '“^ and /i = cos^. Transforming this expression to the emitters’ frame, 
moving at T' = (1 — in the x'-direetion in the jet’s bulk frame (notiee that T' ean be a 

funetion of 7'), we obtain: 

U'^ = (r' 7 '(1 - r' 7 '(/3e/i - /?'), 7e/9e\/r^, 0) (C.2) 


The average momentum in the x' direetion should be zero in the emitters’ frame, as it may be 
defined as the partieles’ loeal eenter of momentum frame. This implies that: 

0 = {u'^) = - 13')) = r'7'(/3e(/i) - f3') P' = m) . (C.3) 


The average (/x) ean be written in terms of 9o or fiQ = cos ~ 1 — l9l'. 


(/^) 


1 

I — fio 



1 + /^o 



(C.4) 


Note that = arccos(/i) 9qI\/2 differs from (9) = {9qcos9q — sin6*o)/(l — cos6*o) 
Given (3' from Eq. |C.3|in terms of the average properties, we ean also obtain T': 


19. 


71 - m? 



ie ieh « 1 , 

^ 7e^0>l. 


(C.5) 


Plugging this baek into the zeroth eoordinate of the four veloeity in the emitters’ frame we get that 
the average energy of partieles in this frame is given by: 

i: = (u«) = r'lKi - m?) = 7. (C.6) 

as given by Eq. 
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